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1 .  INTRODUCTION 


1.1  Opening  Remarks 

From  pre-historic  times  tribological  concepts  and  systems  have 
been  developed  and  used  to  its  advantage  by  mankind.  The  use  of 
lubricants  by  the  Egyptians  in  the  transport  of  large  building  blocks 
and  the  early  attempts  to  develop  rolling-element  bearings  in  Greek 
and  Roman  times  are  good  examples.  It  is,  however,  not  until  the 
mid- twentieth  century  that  many  scientists  and  engineers  began  to 
give  serious  attention  to  the  fundamentals  of  tribology,  although 
the  first  scientific  studies  of  friction  dates  back  to  the  era  of 
Leonardo  da  Vinci  [1]. 

Among  various  theories  proposed  in  the  past,  the  adhesion 
theories  for  friction  and  wear  have  been  widely  accepted.  Most  theories, 
however,  including  the  adhesion  theories,  do  not  predict  some  of 
experimentally  observed  phenomena  and  sometimes  lack  physical  basis. 

In  1973  Suh  [2],  considering  these  problems,  introduced  the  delamination 
theory  of  wear  for  metals  sliding  at  low  speeds.  Extensive  experi¬ 
mental  and  analytical  work  conducted  since  that  time  has  proved  most 
hypotheses  in  the  theory  except  the  exact  mechanism  of  crack  propagation. 
One  of  the  difficulties  in  predicting  the  crack  propagation  rates  was 
the  lack  of  precise  knowledge  on  surface  traction.  The  validity  of 
friction  theories  proposed  in  the  past  has  been  questioned  since  the 
experimentally  observed  high  coefficient  of  friction  could  not  be 
explained  by  the  adhesion  theory  from  the  simple  continuum  mechanics 
point  of  view,  i.e.,  the  yield  condition  is  violated  at  the  interface. 
Therefore,  the  frictional  behavior  of  materials  was  investigated  at 


M.  I.  T.  to  clarify  the  friction  mechanism.  Tohkai  studied  the  effect 
of  varying  hardness  on  friction  in  1979  as  part  of  his  S.  M.  thesis  [3], 
showing  that  the  frictional  behavior  of  materials  generally  affects 
the  wear  behavior,  which  is  consistent  with  the  result  of  the  previous 
investigation  by  Rabinowicz  [4].  This  is  because  of  the  influence  of 
friction  on  crack  propagation  as  well  as  plastic  deformation  and  crack 
nucleation.  Therefore,  it  is  necessary  to  acquire  better  understanding 
of  the  magnitude  of  frictional  force  at  each  asperity  contact  and  then 
to  apply  this  knowledge  in  studying  crack  propagation  so  as  to  determine 
the  wear  rate  of  sliding  metal  surfaces. 

1.2  A  Review  of  Previous  Work  in  Friction  and  Wear 

Since  Amontons  published  the  result  of  his  experimental  investi¬ 
gation  of  the  friction  of  unlubricated  solids  in  1699,  such  investigators 
as  de  la  Hire  (1732),  Euler  (1750),  Coulomb  (1785),  Morin  (1835)  and 
Him  (1854)  contributed  to  the  classic  laws  of  friction  [5].  They 
considered  that  friction  is  due  to  the  interlocking  of  asperities  and 
the  friction  force  is  the  tangential  force  necessary  to  lift  the 
asperities  of  one  surface  over  the  asperities  of  the  other,  which  is 
generally  referred  to  as  the  roughness  theory. 

The  adhesion  theory  of  friction  began  to  attract  the  scientists 
when  Tomlinson  (1929)  and  Hardy  (1936)  attributed  friction  forces  to 
energy  dissipation  when  the  molecules  are  forced  into  each  other's 
atomic  fields  and  then  separated.  At  the  time  the  theory  seemed  to 
explain  the  variation  of  friction  with  varying  surface  contamination 
better  than  the  roughness  theory.  However,  the  adhesion  theory  was 
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established  when  Holm  [6],  Ernst  and  Merchant  [7],  and  Bowden  and 
Tabor  [8]  pointed  out  the  difference  between  the  apparent  area  and  the 
real  area  of  contact,  and  showed  that  the  friction  force  is  independent 
of  the  apparent  area  of  contact,  but  it  is  proportional  to  the  real 
area  of  contact. 

In  early  1940' s,  several  investigators  suggested  quantitative 
models  of  friction.  Ernst  and  Merchant  [7]  have  proposed  that  the 
frictional  force  F  may  be  written  as 

F  =  S  +  R  (1.1) 

where  S  is  the  force  required  to  shear  the  metallic  junctions  and  R 
the  force  required  to  lift  the  asperities  of  one  surface  over  those 
of  the  other.  On  the  other  hand,  Bowden  et  al .  [9]  have  suggested 
that  another  factor  to  be  considered  is  the  force  required  to  displace 
the  softer  metal  from  the  path  of  the  slider,  P,  i.e, 

F  =  S  +  P.  (1.2) 

It  follows  that  S  =  As,  R  =  tane  and  P  =  A'H  where  A  is  the  real  area 
of  contact  of  the  metals,  s  the  shear  strength  of  the  softer  metal, 

0  the  slope  of  asperity  contact.  A'  the  cross-sectional  area  of  the 
grooved  track  and  H  the  hardness  of  the  softer  metal . 

Both  of  the  two  groups  were  concerned  about  the  deviation  of  their 
friction  laws  from  Amontons'  law  and  considered  that,  by  neglecting  R 
and  P  since  they  were  believed  to  be  small  in  comparison  to  S,  the 
coefficient  of  friction  may  be  written  as 
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(1.3) 


which  has  been  the  basic  equation  for  the  dominating  adhesion  theory 
of  friction.  Ernst  and  Merchant  derived  the  expression  for  the  value 
of  shear  stress  to  be  applied  to  a  metal  to  lower  its  melting  point  to 
a  certain  temperature  (usually  room  temperature),  while  Bowden  et  al. 
estimated  s  from  S  =  F  -  P  by  using  their  experimental  results.  How¬ 
ever,  the  hardness-shear  strength  relationship  indicates  that  the 
coefficient  of  friction  in  the  limiting  case  is  approximately  equal 
to  1/6.  This  value  is  much  lower  than  those  measured  in  air  as  well 
as  in  vacuum. 

In  order  to  explain  the  discrepancy  between  the  experimentally 
measured  and  the  theoretically  predicted  friction  coefficients,  the 
adhesion  theory  has  relied  on  the  argument  that  the  area  of  contact 
grows  as  the  tangential  force  is  increased  to  the  point  at  which  gross 
sliding  occurs  [10-12].  When  a  tangential  force  F  is  applied  in  addition 
to  the  normal  load  W  there  is  junction  growth  and  the  area  of  contact 
increases  to  A.  McFarlane  and  Tabor  [10]  assumed  that  the  new  pressure 
p  =  W/A  due  to  the  increase  in  the  area  of  contact  and  the  tangential 
stress  s  =  F/A  satisfy  the  criterion  of  plastic  flow  over  the  contact 
region  of  the  following  form 


n2  4.  c2  _  „2 

P  +  as  =  p^ 


(1.4) 


where  a  is  an  appropriate  constant  (3  for  a  two-dimensional  case)  and 
p^  is  the  normal  pressure  when  s  =  0,  i.e.,  p^^  is  equal  to  the  hard¬ 
ness  H.  If  the  relation  between  p  and  s  in  Eq.  1.4  is  true,  then  it 
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can  be  argued  that  s  can  reach  when  p  -►  0.  Using  this  argument 

Tabor  [12]  explained  how  the  coefficient  of  friction  can  reach  infinity 
theoretically.  However,  the  yield  criterion  in  Eq,  1.4  is  not  correct 
at  all  since  other  stress  components  are  completely  ignored,  for 
example,  component  in  two-dimensional  plane  strain  case  which  is 
necessary  to  satisfy  the  equilibrium  condition.  The  slip-line  field 
solution  by  Green  [13]  indicates  that  the  minimum  value  of  p/H  (or  p/p^) 
is  0.5  for  a  flat  contact.  Also,  Gupta  and  Cook  [14]  have  experimentally 
shown  that  the  average  normal  contact  pressure  for  many  practical 
problems  may  be  obtained  by  the  relation  p/H  =0.5. 

The  adhesion  theory  of  friction  became  theoretically  strong  when 
Green  [15]  analyzed  plastic  junctions  using  slip-line  field  theory  under 
combined  normal  and  shear  stress  for  sliding  between  unlubricated  metals. 
Further,  the  elaboration  of  Bowden  and  Tabor's  friction  theory  using 
the  theoretical  solutions  and  plasticine  models  [13]  has  shown  that  the 
distinctive  feature  of  steady  sliding  is  that,  on  the  average,  the 
surfaces  move  parallel  to  each  other,  and  that  theoretical  plastic 
solutions  for  junctions  provide  qualitative  estimates  for  u  by  con¬ 
sidering  the  life  cycle  of  a  junction.  Also,  strong  junctions  that 
are  sufficiently  ductile  to  survive  until  the  normal  stress  becomes 
tensile  were  regarded  as  a  possible  mechanism  for  a  high  value  of  p 
in  vacuum.  However,  the  insufficient  accuracy  in  the  normal  and 
tangential  forces  exerted  through  a  junction  during  its  life  cycle 
makes  it  impossible  to  estimate  p  quantitatively. 

Gupta  and  Cook  [14,  16]  have  extended  Green's  model  and  developed 
an  analytical  model,  using  a  statistical  distribution  of  spherical 


asperities,  to  study  the  contact  problem.  They  included  the  possibility 
of  some  plastic  junctions  supporting  a  tensile  normal  stress  following 
Green's  suggestion,  and  showed  that  the  effect  of  adhesion  is  unimport¬ 
ant  if  the  tensile  stress  is  less  than  (2  -  3)k,  where  k  is  the  shear 
yield  stress.  The  model  predicts  that  the  number  of  asperity  junctions 
increases  by  the  0.91  power  of  normal  load,  while  the  area  of  each 
individual  junction  grows  only  by  the  0.09  power  with  increase  of  load. 
With  the  statistical  description  of  surface  topography  the  model  was 
found  to  agree  well  with  experimental  results  of  friction  coefficients. 
However,  it  is  the  steady  state  topography  that  is  more  relevant  in 
describing  the  contact  behavior  and  it  is  known  that  this  steady  state 
topography  is  independent  of  the  initial  topography  [17]. 

Recently  some  plasticity  solutions  by  modifying  Green's  solution 
have  been  published  [18,  19].  Among  them  Challen  and  Oxley's  approach 
is  noteworthy:  they  showed  in  their  rubbing  model  how  hard  asperities 
deform  the  surface  of  a  softer  metal  in  the  presence  of  interfacial 
adhesion  [19].  This  may  provide  a  theoretical  basis  for  the  adhesion 
component  of  friction. 

Since  the  adhesion  theory  of  friction  emphasizes  the  importance 
of  adhesion  between  asperities,  a  great  deal  of  attention  in  the  past 
has  been  devoted  to  the  role  of  surface  energy  and  the  mutual  solubility 
of  the  contact  materials.  Rabinowicz  [20,  21]  modified  the  adhesion 
theory  of  friction  by  including  the  effect  of  surface  energy  since  the 
original  adhesion  theory  fails  to  explain  the  large  variation  of 
friction  coefficient.  According  to  the  theory,  high  friction  coeffi¬ 
cients  are  expected  for  sliding  materials  when  the  ratios  of  the  surface 
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energy  to  hardness  are  high.  Further,  he  has  shown  that  metals  with 
a  high  tendency  towards  solid  solution  exhibit  higher  friction  coef¬ 
ficients  than  metals  that  do  not  form  solid  solutions  [22].  Although 
many  experimental  results  were  presented  to  support  the  theory,  some 
of  the  interpretations  are  now  suspect  because  surfaces  are  easily 
contaminated  by  chemisorption  and  physisorption  and  the  chemical  compo¬ 
sition  of  the  surface  is  different  from  that  of  the  bulk  [23]. 

As  can  be  noticed  from  the  review  above,  the  adhesion  theory  of 
friction  has  dominated  other  theories  such  as  the  roughness  and  plowing 
theory,  even  though  an  attempt  was  made  to  combine  the  shearing, 
plowing  and  roughness  components  by  Shaw  and  Macks  [24].  The  early 
roughness  theory  may  be  abandoned,  since  the  raising  and  lowering  of  the 
asperities  of  one  surface  over  those  of  the  other  surface  involves  no 
energy  dissipation.  Friction  is  certainly  a  dissipation  mechanism. 

Early  investigators  were  concerned  about  the  energy  dissipation  during 
sliding  friction.  It  is,  however,  only  in  recent  days  that  the  import¬ 
ance  of  plastic  deformation  has  been  recognized  in  the  study  of  friction 
and  wear. 

The  work  done  by  Liu  [25],  Suh  and  Sridharan  [26]  and  Tsuya  [27] 
prompted  Rigney  and  Mirth  [28]  to  propose  the  following  model  for 
sliding  and  friction  based  on  the  energy  balance  during  plastic 
deformation 

=  wtTe 

L  (1.5) 

where  w  is  the  width  of  the  highly  deformed  region  (approximately  equal 
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to  the  wear  track  width),  t  the  thickness  of  the  region,  t  the  shear 
stress,  e  the  average  strain  per  cycle  and  L  the  applied  load. 

However,  this  model  does  not  provide  any  insignt  into  the  real  cause 
of  friction  and  moreover,  the  value  of  e  is  the  most  difficult  to 
obtain  since  e  is  not  the  steady  state  strain  itself  but  the  strain 
increment  per  cycle. 

From  dislocation  theory  it  is  known  that  plastic  deformation  of 
crystalline  materials  deforming  in  glide  occurs  if  the  momentary 
critical  resolved  shear  stress  is  exceeded  on  any  crystallographic  slip 
system.  Based  on  this  Kuhlmann-Wilsdorf  [29]  recently  proposed  the 
following  expression  for  friction: 

U  «  {1-6) 

where  m^^  is  the  Schmid  factor.  This  expression  has  been  found  to  pre¬ 
dict  a  number  of  effects  which  appear  to  be  consistent  with  observations, 
dependent  on  surface  texture,  surface  hardness,  and  surface  temperature. 
However,  the  estimate  for  w  given  by  Eq.  1.6  may  be  regarded  as  the 
adhesion  component  as  Kuhlmann-Wilsdorf  suggested.  Further,  many 
simplified  assumptions  made  for  the  estimation  of  m^^  for  the  case  of 
polycrystals  remain  to  be  refined  in  regard  of  multiple  slip. 

Although  many  investigations  have  been  conducted  about  friction, 
many  questions  still  remain  unanstrared;  no  extensive  effort  has  been 
made  to  isolate  the  relative  contribution  of  each  mechanism  responsible 
for  friction;  for  example,  the  plowing  component  of  friction  has  been 
assumed  to  be  small  and  ignored  in  the  shade  of  the  dominant  adhesion 
theory.  The  fact  that  the  frictional  behavior  does  not  depend  solely 
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on  adhesion  is  supported  by  the  observation  that  the  frictional  force 
depends  on  the  history  of  sliding.  The  work  done  at  M.  I.  T,  by 
Tohkai  [3]  indicates  that  the  frictional  force  undergoes  significant 
changes  during  the  early  stages  of  sliding  before  reaching  steady  state. 
The  time  dependent  nature  of  frictional  behavior  has  been  largely 
ignored  in  the  past,  although  it  is  a  rich  source  of  information  in 
understanding  friction.  The  difference  between  the  static  and  kinetic 
coefficients  of  friction  has  been  quoted  in  the  literature,  but  the 
time  dependent  nature  of  the  kinetic  coefficient  of  friction  has  not 
been  widely  reported  [30]. 

In  the  field  of  wear  the  adhesion  theory  also  has  been  widely 
accepted.  The  theory  postulates  that  wear  particles  are  formed  when 
the  strong  metallic  junctions  get  sheared  and  become  loose.  Holm  [31], 
Burwell  and  Strang  [32]  and  Archard  [33]  contributed  to  the  theory  in 
formulating  the  adhesive  wear  law.  It  was  the  only  theory  in  sliding 
wear  before  the  delamination  theory  of  wear  emerged  in  1973. 

1.3  The  Delamination  Theory  of  Wear  and  Recent  Developments 

in  the  Theory 

The  delamination  theory  of  wear  describes  the  following  sequential 
(or  independent  if  there  are  pre-existing  subsurface  cracks)  events 
which  lead  to  loose  wear  shear  formation  [34]. 

(1)  When  two  surfaces  come  into  contact,  normal  and  tangential 
loads  are  transmitted  through  the  contact  points.  Asperities  of  the 
softer  surface  are  easily  deformed  and  fractured  by  the  repeated 
loading  action,  forming  small  wear  particles.  Hard  asperities  are  also 
removed  but  at  slower  rates.  A  relatively  smooth  surface  is  initially 


generated,  either  when  these  asperities  are  deformed  or  when  they  are 
removed. 

(2)  The  surface  traction  exerted  by  the  harder  asperities  at  the 
contact  points  induces  incremental  plastic  deformation  per  cycle  of 
loading,  accumulating  with  repeated  loading.  The  increment  of  perma- 
ent  deformation  remaining  after  given  cyclic  loading  is  small  compared 
with  the  total  plastic  deformation  that  occurs  in  that  cycle  due  to 
the  reversal  of  shear  deformation. 

(3)  As  the  subsurface  deformation  continues,  cracks  are  nucleated 
below  the  surface.  Crack  nucleation  very  near  the  surface  cannot  occur 
due  to  the  triaxial  state  of  compressive  loading  which  exists  just 
below  the  contact  region. 

(4)  Once  cracks  are  present  (either  by  crack  nucleation  or  from 
pre-existing  voids  and  cracks),  further  loading  and  deformation  causes 
the  cracks  to  extend  and  propagate,  eventually  joining  with  neighboring 
cracks.  The  cracks  tend  to  propagate  parallel  to  the  surface  at  a 
depth  governed  by  material  properties  and  the  state  of  loading.  When 
the  cracks  cannot  propagate  because  of  the  small  stress  concentration 
at  the  crack  tip  due  to  an  extremely  small  surface  traction  at  the 
asperity  contact,  crack  nucleation  is  the  rate  controlling  mechanism. 

(5)  When  the  cracks  finally  shear  to  the  surface,  long  and  thin 
wear  sheets  delaminate.  The  thickness  of  the  wear  sheet  is  controlled 
by  the  location  of  subsurface  crack  growth,  which  is  controlled  by  the 
normal  and  tangential  loads  at  the  surface. 

A  series  of  experimental  study  conducted  at  M.  I.  T.  [35]  has 
substantiated  the  theory,  showing  the  delamination  process  initiates 
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when  the  subsurface  plastic  deformation  causes  the  nucleation  of  voids 
(Fig.  1.1).  With  further  deformation  these  voids  elongate  and  link  up 
to  form  long  cracks  in  a  direction  nearly  parallel  to  the  wear 
surface  (Fig.  1.2).  At  a  critical  length,  these  cracks  shear  to  the 
surface,  yielding  a  wear  particle  in  the  form  of  a  long  thin  sheet  as 
shown  in  Fig.  1.3.  The  top  surface  of  the  wear  sheet  is  generally 
smooth,  while  the  fractured  surface  is  rough,  often  showing  dimples 
(Fig.  1.4). 

The  mechanism  of  crack  nucleation  in  materials  with  inclusions  or 
hard  second  phase  particles  is  well  known,  although  the  exact  mechanism 
in  single  phase  materials  is  not  clear.  Jahanmir  and  Suh  [36] 
investigated  the  location  of  crack  nucleation,  based  on  the  stress 
criterion,  for  circular  inclusions.  It  indicates  that  voids  can  only 
nucleate  in  a  small  region  below  the  sliding  contact  and  the  depth  of 
this  region  increases  with  increase  of  both  the  normal  load  and  the 
friction  coefficient.  It  further  shows  that  crack  nucleation  takes 
place  readily  in  two-phase  materials,  leading  to  a  conclusion  that  crack 
propagation  controls  the  wear  rate  in  many  engineering  materials. 

The  analysis  of  the  mechanics  of  crack  propagation  using  linear 
elastic  fracture  mechanics  has  shown  that  the  stress  intensity  factors 
reach  maximum  at  a  distance  below  the  surface,  indicating  that  the 
cracks  at  this  depth  propagate  faster  than  others  [37].  Also,  the 
change  in  stress  intensity  factor,  and  therefore  the  crack  propagation 
rate,  increases  with  increasing  coefficient  of  friction.  The  stress 
intensity  factors  calculated,  however,  were  much  smaller  than  the 
threshold  values,  implying  that  the  crack  growth  process  does  not 
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Fig.  1.2  Subsurface  deformation  and  crack  formation  in  iron  solid 
solution  [Ref.  35]. 
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take  place  by  large  macroscopic  crack  propagation  but  rather  by  void 
nucleation  and  the  deformation  of  the  voids  generated.  One  of  the 
reasons  for  the  small  values  is  that  the  method  applied  to  calculate 
the  stress  intensity  factors,  namely,  the  pseudo-stress  field  and  the 
weighting  function  method,  may  be  too  crude  for  such  a  geometrically 
complicated  one  in  sliding  wear. 

Recently,  Hills  and  Ashelby  [38],  Rosenfield  [39],  and  Keer 
et  al.  [40]  have  solved  the  same  subsurface  crack  problem  for  the  stress 
intensity  factors.  The  first  two  papers  have  considered  the  Mode  II 
crack  propagation  and  found  that  the  maximum  stress  intensity  factor 
exists  at  the  compression  region  ahead  of  the  asperity.  On  the  other 
hand,  the  other  two  by  Fleming  and  Suh  and  Keer  et  al.  have  concentrated 
on  the  tensile  region.  Nevertheless,  two  points  are  consistent  with 
each  other:  shear  stresses  dominate  over  normal  stresses  and  the  stress 
intensity  factors  are  too  low  to  cause  crack  growth. 

1.4  Scope  of  Research 

As  reviewed  in  the  previous  section  extensive  experimental  investi¬ 
gation  by  Jahanmir  [41]  has  confirmed  that  wear  sheets  are  indeed 
formed  as  a  result  of  plastic  deformation,  crack  nucleation,  and  crack 
propagation.  As  part  of  his  doctoral  thesis  work  at  M.  I.  T.  he  has 
shown  analytically  that  cracks  can  readily  nucleate  in  the  subsurface. 

In  addition,  the  results  show  that  plastic  deformation  and  crack 
nucleation  are  greatly  affected  by  frictional  force  at  the  asperity 
contact.  The  analyses  of  crack  propagation  by  several  investigators 
have  also  suggested  that  the  higher  the  coefficient  of  friction  is 
the  faster  cracks  propagate. 


It  is  well  known  that  the  coefficient  of  friction  has  a  signifi¬ 
cant  effect  on  the  wear  rate  [3,  4].  This  is  due  to  the  effect  of  the 
surface  traction  on  plastic  deformation,  crack  nucleation  and  crack 
propagation  rates.  Therefore,  it  is  required  that  we  properly  understand 
the  genesis  of  friction  because  the  magnitude  of  the  frictional  force 
at  each  contact  point,  which  will  in  turn  affect  the  crack  nucleation 
and  propagation  rates,  cannot  be  determined  without  this  understanding. 

In  order  to  determine  the  mechanism  of  crack  propagation  a  more 
rigorous  method  such  as  the  finite  element  method  is  required  before 
any  conclusion  is  drawn.  The  methods  used  previously  were  either  too 
approximate  or  limited. 

The  objective  of  this  study  is  to  examine  the  genesis  of  friction 
both  experimentally  and  theoretically,  to  identify  the  mechanism  of 
crack  propagation  through  a  rigorous  analysis  based  on  better  under¬ 
standing  of  the  magnitude  of  surface  traction,  and  to  discuss  the 
implications  of  the  mechanism  of  friction  and  crack  propagation  on  wear. 

A  different  theory  is  proposed  to  explain  the  frictional  behavior 
of  materials  based  on  the  experimental  results  in  Chapter  2.  Each 
mechanism  of  asperity  deformation,  plowing  and  adhesion  is  analyzed 
by  using  plasticity  theory.  In  Chapters  3  and  4,  the  mechanism  of 
crack  propagation  is  discussed  through  the  finite  element  analysis. 

It  first  examines  the  validity  of  the  linear  elastic  fracture  mechanics 
approach  and  then  considers  the  plasticity  aspect  of  subsurface  crack 
growth.  Friction  and  crack  propagation  are  related  to  the  wear  rate  of 
materials  in  Chapter  5.  Finally,  based  on  this  investigation, 
conclusions  are  given  in  Chapter  6. 
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2.  THE  GENESIS  OF  FRICTION 


2.1  Introduction 

Recent  research  has  shown  that  the  frictional  behavior  of  materials 
generally  affects  the  wear  behavior.  This  is  due  to  the  fact  that  high 
friction  not  only  expedites  subsurface  deformation  and  crack  nucleation 
but  also  accelerates  crack  propagation.  Thus  the  exact  load  at  each 
asperity  contact  must  be  known  first  to  assess  its  effect  on  crack  pro¬ 
pagation  and  to  compute  the  crack  propagation  rate.  To  be  able  to  do 
this  the  exact  mechanisms  responsible  for  friction  should  be  identified 
first  of  all . 

In  this  chapter  a  different  theory  of  friction  is  proposed  based 
on  the  experimental  and  theoretical  results  and  quantitative  models 
for  the  friction-generating  mechanisms  such  as  asperity  deformation, 
plowing  by  wear  particles  and  hard  asperities,  and  adhesion  are 
presented. 

2.2  Experiments 

A  series  of  experiments  were  conducted  at  M.  I.  T.  to  study  the 
friction  and  wear  behavior  of  various  combinations  of  the  following 
materials:  Armco  iron,  AISI  1020,  1045,  and  1095  steel.  These  iron 
based  metals  with  differing  carbon  contents  have  large  differences  in 
hard-phase  concentrations  and  hardness.  These  materials  were  chosen 
to  minimize  chemical  differences. 


Armco  iron  was  recrystallized  at  973  K  for  one  hour.  AISI  1020 
1045  and  1095  steels  were  austenitized  at  1173  K  for  15  minutes,  oil 
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quenched,  and  then  tempered  at  673  K  for  one  hour  to  obtain  a  sphero- 
dized  microstructure.  The  hardness  and  the  volume  fraction  of  cementite 
are  listed  in  Table  2.1. 

The  experimental  study  was  initiated  by  Tohkai  [3]  as  part  of  his 
S.  M.  thesis  work  and  he  conducted  half  of  the  tests  reported  in 
Tables  2.2  and  2.3,  that  comprises  the  triangular  part  of  the  matrix 
below  the  diagonal  from  upper  left  to  lower  right.  Later  it  was  extended 
to  this  study. 

The  heat-treated  samples  were  polished  with  4/0  abrasive  paper, 
cleaned  with  trichloroethylene  and  weighed  to  an  accuracy  of  0.01  mg 
before  testing.  Some  OFHC  copper  specimens  were  also  tested  in  air 
against  AISI  1020  steel  to  compare  the  initial  frictional  behavior. 

The  initial  coefficient  of  friction,  was  about  the  same  as  those 
obtained  with  steel  and  iron  specimens.  This  initial  frictional  behavior 
is  being  studied  by  Willett  extensively  [42]. 

Samples  of  6.35  mn  in  diameter  were  tested  for  friction  and  wear 
using  cross-cylinder  geometry.  The  specimen  (rotating  cylinder)  was 
rotated  by  the  spindle  of  the  lathe,  and  the  slider  (stationary  cylinder) 
was  held  stationary  in  a  holder  attached  to  a  lathe  tool  dynamometer 
which  was  mounted  on  the  carriage  of  the  lathe.  Both  normal  and  tan¬ 
gential  forces  were  measured  by  a  dynamometer-recorder  assembly. 

Tests  were  conducted  in  a  purified  argon  atmosphere  except  for 
AISI  1020  steel,  where  some  samples  were  tested  in  air  under  both 
lubricated  and  unlubricated  conditions.  Water  and  light  machine  oil 
were  used  as  lubricants.  The  experimental  results  were  obtained  under 
the  following  conditions:  a  normal  load  of  1  kg  (9.8  N);  a  sliding 
speed  of  0.02  m/s;  a  total  sliding  distance  of  36  m;  room  temperature. 
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Table  2.1.  Experimental  Materials 


Material 

Heat  Treatment 

Vickers 

Hardness 

(MPa) 

Volume 
Fraction  of 
Cementite 

Armco  Iron 

973  K,  1  hr;  air-cooled 

980  +  50 

0.0004 

AISI  1020 
steel 

Austenitized  at  1173  K, 

15  min;  oil -quenched; 

673  K,  1  hr;  air-cooled 

1710  +  100 

0.020 

AISI  1045 
steel 

Spheroidized: 

1173  K,  15  min;  oil-quenched; 
673K,  1  hr;  air-cooled 

4120  +  130 

0.067 

AISI  1095 
steel 

Spheroidized: 

1173  K,  15  min;  oil- 
quenched  >673  K,  1  hr;  air¬ 
cooled 

6080  +  350 

0.142 

Slider  (stationary  cylinder) 
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Table  2.2.  Friction  Coefficients 


specimen  (rotating  cylinder) 

Armco 

1045 

BHjH 

iron 

Hi 

steel 

Armco 

m 

0.13 

0.20 

0.24 

0.20 

iron 

Q 

0.71 

0.75 

0.69 

0.76 

u 

■o 

c 

1020 

■ 

0.18 

0.20 

0.13 

0.12 

>» 

U 

steel 

H 

0.55 

0.68 

0.57 

0.65 

fO 

c 

0.80 

____ 

— - 

o 

•r^ 

4-» 

4-» 

1045 

B 

0.16 

0.17 

0.17 

0.12 

</) 

U 

steel 

B 

0.52 

0.53 

0.71 

0.69 

•o 

0.77 

0.71 

____ 

____ 

1095 

^i 

0.17 

0.17 

0.14 

0.17 

steel 

‘^s 

0.51 

0.54 

0.58 

0.67 

u* 

0.76 

0.73 

=  initial  coefficient  friction 
Ug  =  steady  state  coefficient  of  friction 
M*  =  peak  value  of  the  friction  coefficient 


Specimen  (rotating  cylinder) 
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Some  of  the  specimens  were  sectioned  along  the  sliding  direction 
to  measure  the  slope  of  the  asperity  from  micrographs  since  the  asperi¬ 
ties  of  the  machined  surfaces  were  orientation  dependent. 

In  the  study  of  the  effects  of  ion  implantation  on  friction  and 
wear  by  Shepard  [43],  extremely  well  polished  surfaces  were  slid  against 
each  other  to  investigate  the  surface  damage  after  predetermined  amounts 
of  sliding.  These  surfaces  were  observed  using  scanning  electron 
microscopy. 

2.3  Experimental  Results  and  Observations 

The  friction  and  wear  coefficients  of  the  iron-carbon  system  are 
tabulated  in  Tables  2.2  and  2.3,  respectively.  There  are  several 
important  results  worth  considering  in  detail.  First  of  all,  the 
coefficient  of  friction  changes  as  a  function  of  the  distance  slid, 
especially  at  the  early  stages  of  sliding.  It  usually  has  a  low  initial 
value  and  gradually  increases  until  reaching  a  steady  state  value. 

After  it  reaches  a  maximum  value  the  friction  coefficient  sometimes 
drops  down  if  the  stationary  slider  is  much  harder  than  the  moving 
specimen.  The  same  pair  of  materials  does  not  show  the  drop  in  the 
coefficient  of  friction  when  their  roles  are  reversed.  The  initial 
coefficient  of  friction  is  always  in  the  range  of  about  0.1  to  0.2 
regardless  of  the  materials  tested  and  whether  or  not  lubricants  are 
used.  Second,  the  steady  state  coefficient  of  friction  and  the  wear 
rates  are  higher  when  identical  metals  are  slid  against  each  other  than 
when  a  harder  stationary  slider  is  slid  against  a  moving  softer  specimen. 
However,  when  a  softer  stationary  slider  is  slid  against  a  harder  moving 
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specimen,  the  steady  state  coefficient  of  friction  is  nearly  the  same 
as  those  of  the  identical  materials  sliding  against  each  other.  In 
this  case  the  wear  rates  of  unidentical  pairs  of  materials  are  much 
greater  than  those  of  identical  metals. 

These  changes  in  the  friction  and  wear  behavior  are  related  to  the 
changes  in  the  surface  topography  as  shown  in  Figs.  2.1  and  2.2,  which 
are  the  micrographs  of  the  slider  surface  and  the  specimen  surface, 
respectively.  These  figures  show  that  when  the  stationary  slider  is 
harder  than  the  specimens,  the  hard  surface  is  polished  to  a  mirror 
finish  and  the  high  spots  of  the  softer  surface  also  acquire  the  same 
mirror  finish.  When  the  delamination  process  is  completed  underneath 
the  polished  surface,  new  high  spots  will  be  polished  to  a  mirror  finish. 
This  does  not  happen  when  the  stationary  slider  is  softer  than  the 
specimen  or  when  identical  metals  are  slid  against  each  other.  In  these 
cases  many  plowing  grooves  are  observed  [3]  and  the  surface  always 
stays  rough. 

From  these  experimental  results  the  following  observations  are 

made: 

(1)  The  coefficients  of  friction  vs.  sliding  distance  (or  time) 
may  be  summarized  by  using  two  typical  plots  shown  in  Fig.  2.3.  The 
behavior  of  Fig.  2.3a  always  holds  when  identical  metals  are  slid 
against  each  other.  The  drop  in  the  coefficient  of  friction  in  Fig.  2.3b 
is  associated  with  mutual  polishing  of .the  mating  surfaces  [3,  44]. 

The  behavior  shown  in  Fig.  2.3b  results'  primarily  when  the  hardness  of 

the  stationary  slider  is  much  greater  than  the  moving  specimen. 


Fig.  2.1 
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Scanning  electron  micrographs  of  the  surfaces  of  worn  sliders: 
(a),  (b),  (c),  (d)  iron  on  iron,  1020,  1045,  and  1095  steel, 
respectively;  (e),  (f),  (q),  (h)  1020  on  iron,  1020,  1045,  and 
1095  steel,  respevtively;  (i),  (j),  (k),  (1)  1045  on  iron,  1020, 
1045,  and  1095  steel,  respectively;  (m),  (n),  (o),  (p)  1095  on 
iron,  1020,  1045,  and  1095  steel,  respectively. 
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Fig.  2.2  Scanning  electron  micrographs  of  the  surfaces  of  worn  specimens: 
(a),  (e),  (i),  (m)  iron  on  iron,  1020,  1045,  and  1095  steel, 
respectively;  (b),  (f),  (j),  (n)  1020  on  iron,  1020,  1045,  and 
1095  steel,  respectively;  (c),  (g),  (k),  (o)  1045  on  iron,  1020, 
1045,  and  1095  steel,  respectively;  (d),  (h),  (1),  (p)  1095  on 
iron,  1020,  1045,  and  1095  steel,  respectively. 
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Fig.  2 


.3  Coefficient  of  friction  vs.  sliding  distance:  (a)  for  an  Armco 
iron  slider  sliding  against  an  Armco  iron  specimen  (y.  =  0.13, 
y  =  0.71);  (b)  for  an  AISI  1095  steel  slider  sliding  against 
afi  Armco  iron  specimen  (y^  =  0.17,  y^  =  0.51,  Ayp  =  0.25). 
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(2)  When  wear  particles  are  brushed  from  the  sliding  interface, 
the  coefficient  of  friction  decreases  to  a  low  value  and  gradually 
reaches  a  steady  state  value  again,  as  schematically  illustrated  in 
Fig.  2.4.  (The  effect  of  wear  particles  on  the  friction  coefficient 
was  also  reported  by  Abrahamson  et  al.  [30]  and  more  recently  by 
Kuwahara  and  Masumoto  [45]. 

(3)  The  coefficient  of  friction  can  differ  by  as  much  as  0.2 
even  for  the  same  pair  of  materials  (which  are  chemically  identical) 
depending  on  which  is  a  stationary  slider  and  which  is  a  moving 
specimen  (see  Table  2.2). 

(4)  The  initial  value  of  the  kinetic  coefficient  of  friction 

is  in  the  neighborhood  of  0.1  to  0.2  (but  largely  in  the  range  of  0.12 
to  0.17)  for  many  materials  tested,  i.e.,  gold  on  gold,  steel  on  steel, 
brass  on  steel,  etc.  [3,  44],  and  also  regardless  of  whether  or  not 
lubricants  are  used. 

(5)  When  the  friction  test  is  done  with  extremely  well  polished 
surfaces,  plowing  grooves  are  formed  from  the  onset  of  testing  (see 
Figs.  2.5  and  2.6). 

2.4  Postulated  Genesis  of  Friction 

The  experimental  results  clearly  indicate  that  the  observed 
friction  coefficients  cannot  be  explained  in  terms  of  the  adhesion 
theory  alone.  The  effect  of  entrapped  wear  particles  and  the  existence 
of  which  is  independent  of  environmental  conditions  and  materials 
tested  cannot  be  explained  by  the  adhesion  theory.  The  theory  is 
further  defied  by  the  dramatic  changes  in  the  coefficient  of  friction 
when  the  role  of  the  slider  and  the  specimen  is  reversed. 
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Fig.  2 


.4  Effect  of  removing  wear  particles  for  an  Armco  iron  slider 
sliding  against  an  Armco  iron  specimen  (y.  =  0.13,  u  =  0.71 
y.’  =  0.4).  ’  s 
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Fig.  2.6 


ical  wear  scar  of  an  iron  pin  after  1  revolution  of  sliding 
'icated  test  in  air;  load  =  400  g  [Ref.  43]. 
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Based  on  the  experimental  results  discussed  in  the  preceding 
section  the  following  postulate  is  advanced  to  explain  the  genesis  of 
friction  between  the  sliding  surfaces: 


"The  coefficient  of  friction  between  the  sliding  surfaces 
is  due  to  the  various  combined  effects  of  asperity 
deformation,  plowing  by  wear  particles  and  hard 
surface  asperities,  Up,  and  adhesion  between  the  flat 
surfaces,  p,.  The  relative  contribution  of  these  components 

d 

depends  on  the  condition  of  the  sliding  interface,  which  is 
affected  by  the  history  of  sliding,  the  specific  materials 
used,  the  surface  topography  and  the  environment." 

In  this  section  this  postulate  will  be  described  further. 
Quantitative  analysis  of  each  of  these  friction  components  is  given  in 
later  sections  from  a  theoretical  point  of  view. 

In  order  to  clarify  the  above  postulate,  the  time  dependent  behavior 
of  friction  will  be  considered  qualitatively  by  sub-dividing  it  into 
the  following  stages  (see  Fig.  2.7). 

(a)  Stage  I  --  In  this  early  stage  the  coefficient  of  friction 
seems  to  be  dictated  by  plowing  of  the  surface  by  asperities.  Adhesion 
does  not  play  any  significant  role  in  this  state  due  to  the  contaminated 
nature  of  the  surface.  The  deformation  of  asperities  does  take  place 
at  the  onset  of  sliding  which  affect  the  static  coefficient  of  friction. 
However,  it  appears  that  in  Stage  I,  asperity  deformation  is  not  the 
major  factor  that  determines  the  coefficient  of  friction,  since  they 
deform  as  soon  as  sliding  commences  and  the  surface  is  easily  polished 
with  the  generation  of  new  asperities  in  Stage  I  (see  Fig.  2.8). 
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Deformation  of  rough  asperities  on  1020  steel  specimen:  (a) 
initial  surface;  (b)  after  1  revolution;  (c)  surface  appearance 


Consequently  the  coefficient  of  friction  in  this  state,  ,  is  largely 
independent  of  material  combinations,  the  surface  conditions  and  the 
environmental  conditions. 

(b)  Stage  II  --  In  this  second  stage,  the  frictional  force  begins 
to  rise  slowly  due  to  increase  in  adhesion.  When  the  interface  is 
lubricated.  Stage  I  persists  for  a  long  time  and  Stage  II  may  not  be 
present.  The  slope  in  Stage  II  can  be  steeper  if  the  wear  particles 
generated  by  the  asperity  deformation  and  fracture  are  entrapped  between 
the  sliding  surfaces  and  plow  the  surfaces. 

(c)  Stage  III  -  This  stage  is  characterized  by  a  steep  increase 
in  slope  due  to  the  rapid  increase  in  the  number  of  wear  particles 
entrapped  between  the  sliding  surfaces  as  a  consequence  of  higher  wear 
rates.  The  slope  can  also  be  affected  by  the  increase  in  adhesion  due 
to  the  increase  in  clean  interfacial  areas.  The  force  required  to 
deform  the  asperities  will  continue  to  contribute  to  the  frictional 
force  in  this  stage  as  long  as  surface  asperities  are  present.  The 
wear  particles  are  generated  when  the  process  of  wear  particle  forma¬ 
tion  by  subsurface  deformation,  crack  nucleation  and  crack  propagation 
postulated  by  the  delamination  theory  of  wear  [2]  is  completed.  Some 
of  the  wear  particles  get  entrapped  between  the  surfaces,  causing 
plowing.  The  plowing  will  be  greater  when  the  wear  particles  are 
entrapped  between  metals  of  nearly  equal  hardness,  because  they  will 
penetrate  into  both  surfaces,  preventing  any  slippage  between  the 
particle  and  the  surface. 


(d)  Stage  IV  -  This  stage  is  reached  when  the  number  of  wear 
particles  entrapped  between  the  interface  remains  constant.  This  occurs 
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when  the  number  of  the  newly  entrapped  particles  equal  the  number  of 
entrapped  particles  leaving  the  interface.  The  adhesion  contribution 
to  friction  also  remains  constant  in  Stage  IV.  The  asperity  deform¬ 
ation  is  not  as  important  as  plowing  since  the  asperities  deform 
readily  and  reach  the  equilibrium  state  in  shape  which  is  independent 
of  the  initial  surface  roughness.  When  two  like  metals  are  slid 
against  each  other  or  when  the  mechanism  responsible  for  Stage  V  does 
not  play  a  significant  role,  the  coefficient  of  friction  in  Stage  IV 
is  the  steady  state  frictional  coefficient  between  the  two  metals. 

(e)  Stage  V  —  In  some  cases,  such  as  when  a  very  hard  stationary 
slider  is  slid  against  a  soft  specimen,  the  asperities  of  the  hard 
surface  are  gradually  removed,  creating  a  mirror  finish  as  shown  in 
Fig.  2.9.  In  this  case  the  frictional  force  decreases  due  to  the 
decrease  in  plowing  and  asperity  deformation.  Plowing  decreases  since 
wear  particles  cannot  anchor  on  a  polished  hard  surface. 

(f)  Stage  VI  —  Eventually  when  the  hard  surface  becomes  a  mirror 
smooth  to  a  maximum  extent,  the  softer  surface  also  acquires  the  same 
mirror  finish  and  the  frictional  force  levels  off.  The  surfaces  are 
never  completely  smooth  since  there  are  always  "pot  holes"  due  to  the 
creation  of  delamination  wear  particles.  These  craters  provide  anchoring 
points  for  wear  particles.  When  the  hard  surface  is  not  stationary 

but  moving  against  the  softer  surface,  the  hard  surface  remains  rough, 
probably  because  polishing  cannot  take  place  due  to  geometric  reasons. 

In  this  case.  Stages  V  and  VI  are  not  present. 


2. 5  Analysis  of  the  Friction  Generating  Mechanisms 

The  three  basic  mechanisms  (i.e.,  asperity  deformation,  plowing 
and  adhesion)  that  are  responsible  for  the  generation  of  friction  will 
be  analyzed  in  this  section.  The  asperity  deformation  initially  deter¬ 
mines  the  static  coefficient  of  friction.  Also,  it  affects  the  dynamic 
coefficient  of  friction,  since  asperities  are  continuously  generated 
due  to  delamination  of  wear  sheets.  However,  the  contribution  of  the 
asperity  deformation  to  the  dynamic  coefficient  of  friction  is  not  large 
relative  to  those  by  plowing  and  adhesion,  since  the  asperities  deform 
readily  (perhaps  in  one  asperity  contact),  whereas  the  generation  of 
asperities  has  to  await  the  formation  of  delamination  wear  particles 
which  often  requires  a  large  number  of  cyclic  loading  by  the  asperities 
of  the  opposing  surface.  On  the  other  hand,  plowing  takes  place  con¬ 
tinuously  whenever  wear  particles  are  entrapped  between  the  sliding 
surfaces  or  when  the  asperities  of  the  counterface  plow  in  all  cases 
when  clean  flat  surfaces  come  into  contact  during  steady  state  sliding. 
The  relative  magnitude  of  these  components  will  be  determined  approxi¬ 
mately  by  using  the  slip-line  field. 

2.5.1  Analysis  of  the  Asperity  Deformation, 

Consider  the  two  representative  asperities  approaching 
toward  each  other  as  shown  in  Fig.  2.10.  When  these  asperities  come 
into  contact  with  each  other,  they  have  to  deform  in  such  a  manner 
that  the  resulting  displacement  field  is  compatible  with  the  sliding 
direction  and  that  the  sum  of  the  vertical  components  of  the  surface 
traction  at  the  contacting  asperities  must  be  equal  to  the  applied 
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normal  load.  A  possible  slip-line  field  that  satisfies  the  kinematic 
condition  is  given  in  Fig.  2.11.  The  solution  demands  that  the  shear 
stress  along  OA  be  whatever  is  necessary  to  satisfy  the  condition  that 
a  =  e,  which  is  necessary  to  constrain  the  resulting  deformation  in 
the  sliding  direction.  This  would  be  possible  even  under  the  lubri¬ 
cated  conditions  if  the  interface  OA  is  not  perfectly  smooth  but  rough 
enough  to  allow  mechanical  interlocking.  The  derivation  of  the  normal 
and  tangential  force  corresponding  to  the  slip-line  field  shown  in 
Fig.  2.11  is  given  in  Appendix  A. 

The  general  solution  is  sketched  in  Fig.  2.12.  If  it  is  assumed 
that  the  asperity  deformation  is  the  only  phenomenon  that  takes  place 
at  the  interface  and  is  entirely  responsible  for  the  frictional  force 
Dnder  a  given  load,  the  coefficient  of  friction  due  to  asperity  deform¬ 
ation  Uj  varies  from  0.39  to  1  as  the  slope  of  asperities  increases 
from  0  to  45°.  These  values  are  closer  to  the  static  coefficient  of 
friction  than  the  dynamic  friction  coefficient  measured  during  the 
early  stage  of  sliding,  i.e..  Stage  I  [46]. 

In  dynamic  situations  most  of  the  normal  load  is  carried  by  the 
entrapped  wear  particles  and  the  flat  contacts.  Therefore,  the  actual 
contribution  of  the  asperity  deformation  to  the  frictional  force  is 
expected  to  be  a  small  fraction  of  the  estimated  value.  This  may  be 
estimated  indirectly  from  experimental  results,  as  discussed  later  in 
this  chapter. 

It  should  be  noted  that  the  slip-line  analysis  done  to  determine 
Uj  is  not  very  different  from  those  of  Green  [13],  who  presented 
similar  solutions  to  justify  the  adhesion  theory  by  considering  the 
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deformation  of  asperities  under  the  influence  of  the  load  exerted  by 
welded  asperity  junctions  of  various  strengths.  An  important  difference, 
however,  is  that  in  the  analysis  presented  in  this  chapter  the  asperity 
deformation  is  assumed  to  be  governed  by  the  kinematic  considerations 
regardless  of  the  nature  of  bonding  between  the  contacting  asperities. 
Therefore,  does  not  depend  sensitively  on  adhesion  and  represents 
the  frictional  force  due  to  the  deformation  of  all  interacting  asperi¬ 
ties.  Before  the  onset  of  sliding  between  two  surfaces  determines 
the  static  coefficient  of  friction. 

2.5.2  Analysis  of  Adhesion  Component  of  the  Friction 

Coefficient,  u. 

- ^ - 0^ 

A  frictional  force  can  arise  due  to  the  adhesion  of  two 
nearly  flat  surfaces.  Unlike  the  deformation  of  asperities  this 
frictional  force  is  a  function  of  the  adhesion  between  the  two  opposing 
surfaces.  The  adhesion  force  arises  either  due  to  the  welding  of  two 
nearly  flat  portions  of  the  surface  or  when  the  atoms  are  brought 
together  to  close  proximity  for  interatomic  interactions  but  without 
welding.  The  adhesion  at  the  slopes  of  two  interacting  asperities  may 
be  neglected  in  this  analysis  since  they  have  already  been  considered 
in  deriving  in  Also  they  are  kinematically  constrained  deformations, 
being  independent  of  the  quality  of  adhesion. 

Experimental  results  show  that  is  not  present  (or  is  at  least 

O 

negligible)  at  the  onset  of  sliding,  probably  due  to  the  presence  of 
contaminants  on  the  surface.  With  the  deformation  of  asperities  and 
exposure  of  fresh  new  surfaces  the  adhesion  between  nearly  flat  surfaces 
is  expected  to  increase.  The  exact  adhesion  area  cannot  be  determined 
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a  priori,  since  the  applied  normal  load  may  also  be  carried  by  inter¬ 
acting  asperities  and  entrapped  wear  particles,  although  the  limiting 
cases  can  be  analyzed. 

Consider  two  nearly  flat  surfaces  coming  into  contact  as  shown 
in  Fig.  2.13.  (Sometimes  this  type  of  contact  is  called  a  "rubbing" 
contact).  Depending  on  the  nature  of  adhesion  along  the  interface  ED, 
the  force  required  to  move  the  interacting  surfaces  with  respect  to  each 
other  will  vary.  When  there  is  no  adhesion  the  force  will  be  zero  and 
when  there  is  complete  adhesion  it  will  reach  a  maximum.  The  solution 
to  this  problem  can  be  obtained  again  using  the  slip-line  fields 
similar  to  that  shown  in  Fig.  2.13.  The  exact  geometric  shape  of  the 
slip-line  field  will  depend  on  the  boundary  condition  at  ED.  The 
solution  sought  can  be  adapted  from  the  recent  work  of  Challen  and 
Oxley  [19]  who  derived  an  expression  for  the  friction  coefficient  as 

_  A  sin  a  +cos  (cos“V-a) 

“a  -  ;  ,1 - (2.1) 

A  cos  o  +sin  (cos  f-o) 


where 


A  =  1  +  I-  +  cos"'  f  -  2o  -  2  sin"'  (2.2) 

^  /TT 


f  ■  strength  of  the  adhesion  at  ED  as 
expressed  at  a  fraction  of  the  shear 
flow  of  strength  of  the  softer  material 

a  *>  the  slope  of  the  hard  asperity 
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For  nearly  flat  surfaces  a  0.  Therefore,  Eqs.  2.1  and  2.2  reduce 
to 


f 

- n — 

®  A  +  sin  (cos  f) 

A  =  1  2  cos'^  f 

varies  from  0  to  0.39  as  f  changes  from  0  to  1. 

The  friction  coefficient  determined  by  Eqs.  2.1  and  2.2  are  based 
on  the  assumption  that  all  the  applied  normal  load  is  carried  by  the 
flat  interfaces.  However,  since  part  of  the  normal  load  is  also  carried 
by  purely  elastic  contacts,  the  interacting  asperity  junctions  discussed 
in  the  preceding  section,  and  the  entrapped  wear  particles,  under 

a 

typical  sliding  conditions  should  be  less  than  0.4.  The  experimental 
results  obtained  with  the  hard  AISI  1095  steel  slider  and  the  soft 
Armco  iron  specimen  showed  that  the  steady  state  coefficient  of  friction 
reached  a  value  of  0.51  when  both  surfaces  were  polished  smooth  and 
thus  the  friction  was  caused  primarily  by  adhesion.  The  agreement 
between  the  theory  and  the  experiment  is  reasonable  since  asperity 
interactions  and  plowing  by  wear  particles  must  have  also  contributed 
to  the  frictional  force. 

2.5.3  Analysis  of  Plowing,  Pp 

The  plowing  component  of  the  frictional  force  can  be  due 


(2.3) 

(2.4) 


to  the  penetration  of  hard  asperities  or  due  to  the  penetration  of 
wear  particles.  The  plowing  due  to  wear  particles  is  schematically 
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illustrated  in  Fig.  2.14.  When  two  surfaces  are  of  equal  hardness, 
the  particle  can  penetrate  into  both  surfaces.  As  the  surfaces  move 
with  respect  to  each  other,  grooves  will  be  formed  in  one  or  both  of 
the  surfaces.  When  one  of  the  surfaces  is  very  hard  and  smooth,  the 
wear  particle  will  simply  slide  along  the  hard  surface  and  no  plowing 
can  occur.  However,  when  the  hard  surface  is  very  rough,  wear  particles 
can  anchor  in  the  hard  surface  and  plow  the  soft  surface. 

The  friction  due  to  plowing  was  investigated  by  Sin  et  al.  [47], 
who  showed  that  the  contribution  of  plowing  to  the  friction  coefficient 
is  very  sensitive  to  the  ratio  of  the  radius  of  curvature  of  the  par¬ 
ticle  to  the  depth  of  penetration.  The  plowing  force  may  be  expressed 
as  P  =  pAg,  where  Ag  is  the  cross-sectional  area  of  the  grooved  track 
and  p  the  flow  stress.  Considering  that  the  contact  interface  is  only 
on  the  front  half  of  the  particle  [48],  we  may  derive  the  expression 
for  the  friction  coefficient  by  plowing  Up  as  (Appendix  B): 


Wp  =  f  [(^)'  sin-^ 


—  -  {(^)^  - 
2r  ^^w^ 


(2.5) 


where  w  is  the  width  of  the  penetration  and  r  is  the  radius  of  curvature 
of  the  particle.  The  ratio  of  w/r  measured  by  sectioning  the  worn 
specimen  is  in  the  neighborhood  of  0.8.  Substituting  this  value  into 
Eq.  2.5,  the  plowing  coefficient  of  friction  is  found  to  be  0.2.  This 
value  is  in  the  same  range  as  the  decrease  in  the  friction  coefficient 
observed  by  removing  the  wear  particles  from  the  Armco  iron/Armco  iron 
and  Armco  iron/AlSI  1995  steel  interfaces,  which  were  0.31  and  0.16 
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respectively.  The  range  of  possible  values  of  Up  as  a  function  of 
the  ratio  w/2r  is  shown  in  Fig.  2.15. 

Plowing  not  only  increases  the  total  frictional  force  and  delam¬ 
ination  wear,  but  also  creates  small  wear  particles  which  in  turn 
affect  the  subsequent  wear  of  sliding  surfaces.  Plowing  action  forms 
ridges  along  the  sides  of  plowed  grooves.  When  these  ridges  are 
deformed  flat  and  subjected  to  repeated  loading,  some  of  these  become 
loose  wear  particles  with  continued  sliding.  This  is  schematically 
illustrated  in  Fig.  2.16. 

2.6  Relative  Contributions  of  u.,  u.,,  and  u  to  the  Overall 

Frictional  Behavior  ^  P 

If  the  postulate  for  friction  coefficient  presented  in  this  chapter 
is  correct,  then  the  relative  values  of  various  friction  components 
are  as  follows: 

Ud  “■  0*43  to  0.75  when  the  entire  applied  normal  load  is 

carried  by  typical  surface  asperities  with  a  slope  of  4° 
to  20°.  It  appears  that  is  responsible  for  the  static 
coefficient  of  friction,  but  does  not  contribute  to  the 
initial  coefficient  of  friction  in  Stage  I.  The 
reason  y^  is  not  a  major  factor  in  Stage  I  is  that  once  the 
original  asperities  deform  asperity  interactions  cannot 
take  place.  This  friction  component  can  contribute  partially 
to  the  steady  state  coefficient  of  friction  if  new  asperities 
are  continuously  generated  as  a  consequence  of  delamination 


wear  process. 


Fig.  2.] 5  Plowing  component  of  the  frict 
function  of  thr  ratio  of  the  w 


Groove 


(b) 


Fig.  2.16  Schematic  illustration  of  wear  particle  formation  due 
to  plowing:  (a)  ridges  formed  along  the  sides  of  the 
plowed  groove;  (b)  flattened  ridges. 
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li  --  From  about  0  to  0.4  depending  on  the  nature  of  adhesion 

a 

between  the  flat  part  of  the  interacting  surfaces.  The 
low  value  is  for  a  well  lubricated  surface  with  light 
lubricant,  while  the  high  value  is  for  identical  metals 
sliding  against  each  other  without  any  surface  contaminants 
and  oxide  layers. 

Up  --  From  nearly  0  to  1.0  from  a  theoretical  point  of  view 

depending  on  the  depth  of  penetration,  but  normally  less 
than  0.4  in  a  typical  situation.  The  high  values  are 
associated  with  two  identical  metals  sliding  against  each 
other  with  deep  penetration  by  wear  particles,  while  the  low 
value  is  obtained  when  either  wear  particles  are  totally 
absent  from  the  interface  or  a  soft  surface  is  slid  against 
a  hard  surface  with  a  mirror  finish. 

The  determination  of  the  total  friction  coefficient  in  a  given 

condition  can  be  rather  complex.  It  is  difficult  to  determine  the 

relative  contributions  of  u.,  u  and  to  the  total  frictional 

da  p 

coefficient,  because  analyses  for  and  were  done  assuming  that 

the  total  normal  load  is  carried  by  either  asperities  or  flat  contact 
areas.  In  real  situations  the  normal  load  will  be  apportioned  among 
the  asperity  contacts,  flat  adhesion  junctions,  and  the  entrapped 
particles.  However,  it  is  quite  plausible  in  many  real  situations  that 
each  one  of  these  mechanisms  which  contribute  to  friction  can  take  place 
sequentially  rather  than  concurrently.  Consider,  for  example,  a  flat 
junction  and  an  asperity  in  contact  as  shown  in  Fig.  2.17.  When  the 
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Fig.  2.17  Two  sliding  surfaces  in  contact:  (a)  first  contact 
of  flat  surfaces;  (b)  flat  surface  contact  and 
asperity  contact. 
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flat  areas  come  Into  contact  first  and  form  an  adhesion  junction,  the 
analysis  performed  for  is  strictly  valid.  When  the  asperities  also 

Q 

come  into  contact  with  further  sliding,  a  large  fraction  of  the  normal 
load  is  still  carried  by  the  flat  junctions  due  to  their  higher  normal 
"stiffness,"  i.e.,  the  force  required  to  cause  unit  displacement  along 
the  vertical  direction.  In  this  case  in  order  to  continue  the  sliding 
action  the  asperities  will  simply  have  to  shear  along  the  dotted  line 
if  the  materials  are  identical  or  along  the  crossed  line  if  the  top 
slider  is  much  harder  than  the  bottom  slider.  In  this  case  only  a 
tangential  force  is  required  to  deform  the  asperity  plastically  along 

the  direction  of  sliding. 

The  coefficient  of  friction  y  is  represented  in  a  "Friction 
Space"  as  a  function  of  adhesion,  plowing  and  roughness  as  shown  in 
Fig.  2.18.  The  adhesion  is  expressed  in  terms  of  the  non-dimensional 
interfacial  shear  strength  f  of  the  flat  contacts  which  was  defined 
by  Eq.  2.2.  The  roughness  is  plotted  in  terms  of  the  slope  of  the 
surface,  while  the  plowing  is  given  in  terms  of  the  ratio  of  the  width 
of  wear  particle  (or  hard  asperity)  penetration  to  its  radius.  The 
lowest  surface  corresponds  to  the  the  case  of  no  asperities,  i.e., 

0  =  0,  which  forms  the  lower  bound.  The  0^  surface  corresponds  to  the 
initial  machined  surface,  while  the  0*  surface  represents  the  surface 
roughness  when  the  peak  of  y  occurs.  As  the  surface  gets  rougher  and/or 
the  number  of  the  steady  state  asperities  increases,  y  will  increase 
and  the  friction  value  will  move  in  the  Friction  Space  along  the  y-axis. 
The  friction  surfaces  0*  is  plotted  from  the  actual  experimental 
results  obtained  with  Armco  iron  sliding  against  AISI  1095  steel. 
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Lower  LIm 
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Fig.  2.18 


"Friction  Space"  showing  the  coefficient  of  friction 
as  a  function  of  adhesion  between  flat  contacts,  wear 
particle  penetration,  and  surface  roughness:  f  =  t  /k 
where  t  is  the  shear  stress  at  the  interface  and  ^ 
k  is  the  shear  flow  strength  of  the  soft  metal;  6  is 
the  slope  of  asperities;  w/2r  is  the  ratio  of  the 
width  of  asperity  penetration  to  the  diameter  of  the 
particle. 


In  this  case,  the  asperity  contribution  to  friction  was  0.3 

Although  this  figure  is  not  precise,  it  gives  a  reasonable  picture 
of  what  happens  in  a  given  situation.  The  paths  of  the  friction  coef¬ 
ficient  change  shown  in  Fig.  2.3  are  shown  in  this  Friction  Space  in 
Fig.  2.18.  The  friction  is  shown  to  start  from  the  initial  roughness 
plane  0^  to  nearly  flat  surface,  i.e.,e=0  and  traces  along  the  paths 
indicated  in  Fig.  2.3.  In  many  cases  the  6^  surface  may  be  above  the 
0*  friction  space,  depending  on  the  initial  surface  finish  relative  to 
the  steady  state  surface  roughness. 

The  foregoing  argument  may  be  applied  to  the  specific  case  of 

gold  sliding  against  gold.  When  unlubricated  gold  specimens  without 

any  oxide  layer  are  slid  against  each  other,  the  frictional  force  is 

due  to  the  sum  of  u.  and  u  or  the  sum  of  a  fraction  of  u u  and 
«  P  da 

u_.  depending  on  the  situation.  It  is  most  likely  that  the  w.  term 
will  always  be  present  since  two  flat  surfaces  in  contact  may  be  the 
most  stiff  system  in  supporting  the  normal  load.  For  example,  if  the 
normal  load  is  first  borne  by  Tat  contacts  only,  friction  will  be 
entirely  due  to  adhesion,  which  will  reach  a  maximum  value.  Then 

a 

if  asperities  of  the  opposing  surfaces  come  in  contact,  an  additional 
frictional  force  will  be  required  to  deform  the  asperities.  In  addition 
to  these  frictional  forces  the  third  frictional  component  may  also 
affect  the  frictional  behavior  if  the  wear  particles  become  wedged  in 
between  the  sliding  surfaces.  Therefore,  the  coefficient  of  friction 
between  gold  on  gold  may  be  as  high  as  1.4  to  1.8  and  fluctuate  between 
a  maximum  and  a  minimum  value. 


70 


Lubricated  'urfaces  can  have  a  coefficient  of  friction  whose 
magnitude  will  be  determined  by  the  degree  of  plowing  and  asperity 
interaction.  Lubricated  surfaces  without  any  asperities  are  found  to 
have  a  coefficient  of  friction  of  approximately  0.04  for  a  hard  surface 
and  0.12  for  a  soft  iron  surface.  However,  when  wear  particles  are 
entrapped  between  the  sliding  surfaces,  the  plowing  component  of  the 
frictional  force  can  be  present  raising  the  friction  coefficient. 

In  the  past  the  high  friction  coefficient  between  like  metals  has 
been  explained  in  terms  of  greater  adhesion  due  to  their  greater 
solubility  [49].  However,  the  evidence  presented  in  this  paper  shows 
that  the  so-called  compatibility  of  metals  is  dictated  more  by  their 
mechanical  behavior  rather  than  by  their  chemical  behavior.  This  is 
quite  reasonable  since  the  diffusion  rate  at  the  typical  sliding 
junctions  is  so  low  due  to  its  low  temperature  that  the  solubility 
between  the  metals  cannot  account  for  the  observed  wear  rates  [50]. 

In  their  classical  book  on  tribology,  Shaw  and  Macks  also  con¬ 
sidered  the  possibility  that  the  frictional  behavior  is  affected  by 
the  asperities,  plowing  and  adhesion.  However,  the  details  of  the 
mechanisms  proposed  in  this  paper  and  their  work  differ  substantially. 

In  their  book,  the  asperity  was  assumed  to  increase  the  friction  due 
to  the  vertical  motion  of  asperities  as  they  slide  over  each  other; 
the  plowing  of  the  sliding  surfaces  by  wear  particles  was  not  con¬ 
sidered  and  the  importance  of  the  penetration  depth  was  not  specifically 
included;  the  adhesion,  plowing  and  asperity  contributions  were  assumed 
to  be  always  constant;  and  the  sliding  distance  dependent  phenomenon 
was  not  considered.  The  work  reported  in  this  paper  clearly  shows 
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(or  hypothesizes)  that  the  adhesion  between  sliding  surfaces  Is  not 
always  present;  that  the  deformation  of  asperity,  rather  than  a  simple 
sllding-over  the  asperities  affects  the  frictional  behavior,  which  Is 
kinematically  constrained;  that  wear  particles  affect  the  coefficient 
of  friction  more  significantly  than  adhesion  In  many  cases;  and  that 
the  coefficient  of  friction  Is  not  a  constant  but  rather  moves  about 
In  the  "Friction  Space"  during  sliding. 

2.7  Sumnary 

The  coefficient  of  friction  Is  shown  to  be  composed  of  three 
components:  due  to  the  deforming  asperities;  Hp  due  to  plowing  by 

wear  particles  and  hard  surface  asperities;  due  to  adhesion.  The 

d 

contribution  to  the  overall  coefficient  of  friction  by  plowing  and 
asperity  deformation  can  be  greater  than  that  by  adhesion,  which  has 
been  largely  Ignored  by  the  classical  adhesion  theory.  Thus,  to 
assess  the  effect  of  friction  on  crack  propagation  and  therefore  on 
wear  all  the  mechanisms  should  be  Included  to  determine  the  load  at 
each  asperity  contact. 
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3.  THE  LIMITATIONS  OF  THE  LEFM  APPROACH 
TO  SUBSURFACE  CRACK  PROPAGATION 


3.1  Introduction 

As  reviewed  previously  five  processes  are  involved  in  delamination 
wear:  transmission  of  forces,  deformation,  crack  nucleation,  crack 
propagation,  and  wear  sheet  formation.  The  wear  rate  of  material  is 
generally  determined  by  the  slowest  or  rate  controlling  process  among 
these  steps.  The  recent  work  on  the  delamination  theory  of  wear  has 
suggested  that  in  many  engineering  materials  the  critical  process  is 
crack  propagation  [36,  37,  41]. 

3.2  Crack  Propagation  Studies  in  the  Past 

Once  cracks  have  been  nucleated  in  the  subsurface  of  sliding 
materials,  they  must  propagate  to  produce  wear  sheets.  Without  crack 
propagation  sheet-like  wear  particles  may  not  form.  The  nucleation  and 
growth  of  several  cracks  may  be  involved  in  the  formation  of  a  single 
wear  sheet.  Sometimes  the  growth  and  coalescence  of  a  large  number  of 
neighboring  voids  will  lead  to  the  formation  of  a  wear  particle. 

Fleming  and  Suh  [37]  have  analyzed  the  propagation  of  a  subsurface 
crack  parallel  to  the  surface  using  a  linear  elastic  fracture  mechanics 
(LEFM)  model.  This  treatment  was  based  on  the  fact  that  a  subsurface 
crack  is  closed  in  the  plastic  region  in  front  of  the  asperity  contact 
and  a  portion  of  the  crack  is  subjected  to  tensile  stress  in  the  elastic 
region  behind  the  contact  zone.  The  stress  intensity  factors  were 
calculated  by  an  elliptical  distribution  of  loads  at  the  contact  using 
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an  approximate  method  based  on  weighting  factors.  A  model  was 
constructed  for  the  correlation  of  fatigue  crack  propagation  and  wear 
rates.  However,  the  correlation  Is  poor  since  the  stress  Intensity 
factors  calculated  are  too  low  to  cause  crack  growth. 

Following  the  Fleming  and  Suh's  work  others  such  as  Hills  and 
Ashelby  [38]  and  Rosenfleld  [39]  have  attempted  to  apply  the  fracture 
mechanics  approach  to  wear.  Both  approaches  are  quite  similar  and  they 
used  the  slightly  different  boundary  loading  conditions  from  that  of 
Fleming  and  Suh.  While  Fleming  and  Suh  concentrate  on  the  tensile 
stress  zone  right  behind  the  contact.  Hills  and  Ashelby  and  Rosenfleld 
consider  the  possibility  of  Mode  II  propagation  of  a  crack  In  a  com¬ 
pressive  field. 

In  considering  a  crack  subjected  to  a  combined  shear-compression 
stress  they  have  assigned  a  friction  coefficient  between  the  crack 
surfaces  In  a  manner  suggested  by  Swedlow  [51].  In  addition  to  stress 
Intensity  calculation  at  a  crack  tip  Hills  and  Ashelby  have  tried  to 
explain  a  possible  mode  of  crack  propagation  under  this  condition  fol¬ 
lowing  McClIntock  [52].  On  the  rail  cracking  problem  McClIntock  has 
used  a  boundary  Integral  relaxation  method  to  calculate  the  plasticity 
at  the  tip  of  a  small  horizontal  crack  In  the  rail  head.  A  wheel 
passage  was  shown  to  give  Initial  sliding,  followed  by  locking, 
sqeezing  the  plastic  zone,  reversed  sliding,  locking,  and  finally 
unloading. 

Recently,  Keer  et  al.  [40]  have  applied  the  Integral  equation 
method  to  an  elastic  half-space  containing  a  horizontal  subsurface 
crack  and  a  surface  breaking  vertical  crack.  They  were  loaded  by  the 
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applied  stress  due  to  the  asperity  contact.  Using  these  models  they 
have  assessed  the  tendencies  of  the  cracks  to  propagate,  to  change 
direction,  and  to  interact  with  each  other.  For  a  horizontal  subsurface 
crack  the  cracks  only  in  a  tensile  field  have  been  studied. 

Although  the  methods  and  the  location  of  cracks  in  each  study 
appear  to  be  different,  one  thing  is  common:  the  predicted  values  of 
stress  intensity  factor  are  too  low  to  cause  crack  growth.  Since  sub¬ 
surface  delamination  cracks  appear  to  grow  under  combined  shear- 
compression  loads,  it  may  not  be  relevant  to  relate  to  fatigue  crack 
propagation  rates  in  Mode.  I.  However,  the  experimental  data  on  Mode  II 
and  combined  Modes  I  and  II  crack  growth  are  not  currently  available 
for  most  materials.  Nevertheless,  it  is  unlikely  that  the  loading 
condition  would  significantly  change  the  crack  propagation  behavior. 
Moreover,  a  recent  study  by  Ritchie  [53]  has  shown  that  wear  resistance 
in  an  ultra  high  strength  steel  (300-M)  bears  no  direct  relationship  to 
crack  propagation  resistance,  whereas  a  direct  correlation  is  found  with 
resistance  to  crack  initiation.  Therefore,  it  is  still  premature  to 
relate  subsurface  crack  propagation  to  fatigue  crack  propagation.  At 
least  following  questions  require  further  investigation: 

(1)  The  accuracy  in  the  methods  of  stress  intensity  factor 
calculation. 

(2)  The  appropriateness  of  the  linear  elastic  fracture  mechanics 
(LEFM)  approach  in  terms  of  plasticity  aspect. 

In  this  chapter  the  stress  intensity  factors  will  firstly  be 
calculated  using  the  finite  element  method.  This  will  include  the 
discussions  on  the  criteria  for  crack  growth  under  sliding  and 
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compression  and  on  the  direction  of  extension  of  cracks.  Then,  the 
plasticity  aspect  at  the  crack  tips  will  be  considered  to  check  the 
appropriateness  of  the  LEFM. 

3.3  Stress  Intensity  Factors  by  Finite  Element  Method 

Various  analytical  techniques  have  been  employed  in  recent  years  to 
determine  stress  intensity  factors  for  through  or  edge  cracks,  and  for 
cracks  embedded  in  infinite  or  semi-infinite  elastic  solids.  However, 
the  analytical  solutions  are  not  available  for  most  cracked  surfaces 
except  for  a  limited  number  of  simple  problems.  Therefore,  the  finite 
element  method  is  suggested  as  the  best  candidate  for  obtaining  stress 
intensity  factors,  whenever  exact  solutions  are  not  available. 

The  finite  element  method  is  one  of  the  most  powerful  and  rigorous 
techniques  available  in  computing  stress  fields,  especially  when  the 
load  applied  is  arbitrary  and  the  geometric  shape  is  irregular.  Since 
the  method  is  conceptually  simple,  applicable  to  large  classes  of 
geometries,  materials  and  loading  conditions,  and  can  be  made  quite 
accurate,  finite  element  methods  are  rapidly  adopted  in  structural 
analysis.  The  details  of  the  method  are  described  in  the  literature  [54] 
and  will  not  be  repeated  here. 

3.3.1  Finite  Element  Model 

3. 3. 1.1  Crack  Tip  Elements 

Recently  several  authors  have  employed  the  8-noded 
and  20-noded  isoparmetric  elements  in  near  tip  modeling  of  2-dimensional 
and  3-dimensional  cracks.  Menshell  and  Shaw  [55]  and  Barsoum  [56]  have 
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demonstrated  that  when  the  mid-side  nodes  in  these  elements  are  placed 
at  the  quarter  point,  the  elements  will  have  the  inverse  square-root 
singularity  characteristic  of  linear  fracture  mechanics.  The  latter 
author  used  both  the  quadratic  quarter  point  element  and  the  triangular 
element  formed  by  collapsing  one  side  (Fig.  3.1),  and  showed  that  the 
triangular  element  leads  to  far  better  results. 

Later  it  has  been  shown  by  Barsoum  [57]  that  the  triangular  elements 
possess  the  same  singularity  in  the  interior  as  well  as  on  the  boundary. 

On  the  other  hand,  quadratic  elements  have  the  singularity  only  on  the 
boundary.  In  addition,  when  two  corner  nodes  and  a  midpoint  node  in  an 
8-noded  element  are  collapsed  in  such  a  way  that  they  initially  are  at 
the  same  location,  but  are  allowed  to  spearate  when  the  element  deforms, 
then  this  element  has  a  1/r  or  1/r  +  1//F  singularity  depending  on  whether 
the  side-nodes  are  at  the  midpoints  or  at  the  quarter  points. 

The  computational  procedure  for  stress  intensity  factors  by  the 
finite  element  method  is  described  in  detail  in  Appendix  C.  It  is 
shown  that  the  stress  intensity  factors  can  directly  be  determined  from 
the  stresses  calculated  by  the  FEM  when  the  triangular  isoparametric 
elements  with  the  mid-size  nodes  placed  at  the  quarter  point  position 
are  employed  as  crack  tip  elements.  In  this  method  both  near  tip  stress 
fields  for  the  crack  tip  element  and  for  the  theoretical  stress  function 
are  compared. 


3. 3. 1.2  Finite  Eletr»nt  Meshes 

The  model  used  for  sliding  wear  study  is  that  of  an 


infinite  half-plane,  two-dimensional  solid  in  plane  strain,  loaded  by 
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1/yr  -  singulority 


1/r  -singularity 


Fig.  3.1  Singular  finite  elements. 
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a  stationary,  rigid  and  straight-sided  asperity.  A  cutout  of  the  semi¬ 
infinite  plane  along  ABCD  was  used  (Fig.  3.2).  The  material  and  loading 

constants  used  were  as  follows:  Young's  modulus,  E  =  1.96  x  10^  MPa  = 

4  2 

2  X  10  kg/mm  ;  Poisson's  ratio,  v  =  0.28;  maximum  normal  stress  applied, 

2 

Pjj  =  980  MPa  =  100  kg/mn  .  This  material  corresponds  to  steel. 

Fig.  3.3  shows  the  geometry  of  the  problem  where  a  crack  of  length 
2c  lies  parallel  to  the  sliding  direction  of  the  asperity  contact  of 
length  2a  and  below  the  surface  at  a  depth  of  d.  Three  different  values 
of  crack  length  were  studied:  c  =  l/4a,  a,  3a.  Figs.  3.4  to  3.6  show 
the  finite  element  meshes  for  each  case  containing  154  elements  and 
498  nodes,  146  elements  and  414  nodes,  and  144  elements  and  488  nodes, 
respectively. 

3.3.2  Results  of  Finite  Element  Calculation 

In  order  to  investigate  the  effects  of  asperity  contact  moving 
along  the  surface  on  subsurface  cracking,  a  number  of  parametric  study 
was  conducted.  For  the  cases  of  c  =  l/4a  and  3a  shown  in  Figs.  3.4 
and  3.6,  respectively,  the  depth  of  crack,  d,  and  the  friction  coeffi¬ 
cient,  VI,  at  the  asperity  contact  were  varied  under  the  moving  asperity 
contact. 

When  a  crack  is  located  in  a  tensile  field,  it  is  expected  to  open 
as  described  by  Fleming  and  Suh  [37].  The  displacement  field  calculated 
for  the  case  of  c  =  a  shown  in  Fig.  3.5  shows  that  the  crack  is  really 
in  opening  mode.  If  the  asperity  moves  to  left,  then  part  of  the  crack 
surfaces  will  tend  to  contact  each  other.  When  this  happens  a  special 
treatment  for  crack  faces  is  required.  Without  any  constraint  along 


Fig.  3.2  A  model  of  a  subsurface  crack  in  sliding  contact  used 
for  the  finite  element  analysis. 


Fig.  3.4  Finite  element  mesh  for  c  =  ]/4a. 
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the  crack,  the  upper  surface  will  displace  down  below  the  lower  surface 
when  loading  since  the  material  over  the  crack  is  connected  to  the  whole 
body  only  at  both  ends  of  the  crack.  To  prevent  this,  the  vertical 
displacements  of  the  upper  crack  surface  were  constrained  lest  it  should 
move  down  below  the  lower  face.  Although  we  may  impose  the  friction 
boundary  condition  along  the  crack  surfaces  (see  Appendix  D),  only  the 
displacement  boundary  condition  was  considered  in  the  study. 

In  Figs.  3.7  to  3.18  the  Mode  II  stress  intensity  factors  for 
c  =  l/4a  and  3a  are  plotted  as  a  function  of  crack  tip  location.  For 
a  given  depth  of  crack  Kjj  increases  to  reach  some  maximum  value  and 
then  decreases  as  the  asperity  moves  towards  the  right.  From  the  figures 
the  following  results  can  be  deduced: 

(1)  The  maximum  values  of  Kjj  occur  when  the  crack  tips  lie  right 
below  the  asperity  contact. 

(2)  The  closer  the  crack  is  to  the  surface,  the  larger  become 
the  maximum  values  of  K^j. 

(3)  A  longer  crack  has  larger  values  of  Kjj. 

These  results  cannot  directly  be  compared  with  the  results  of 
Keer  et  al.  [40]  since  they  have  studied  cracks  only  in  a  tensile  region 
for  subsurface  horizontal  cracks.  Nevertheless,  their  results  show  that 
cracks  even  in  the  tensile  region  have  the  same  properties  as  those  of 
(2)  and  (3)  listed  above.  However,  this  study  indicates  that  cracks  in 
the  tensile  region  have  much  smaller  values  of  Kjj. 
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Fig.  3.10 


Normalized  Mode  II  stress  intensity  factor  vs.  distance 
from  asperity  centerline:  c  =  l/4a;  y  =  0.5;  right  tip. 


centerline 


Normalized  Mode  II  stress  intensity  factor  vs.  distance 
from  asperity  centerline:  c  =  3a;  y  =  0.5;  left  tip. 
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Fig.  3.16  Normalized  Mode  II  stress  intensity  factor  vs.  distance 
from  asperity  centerline:  c  =  3a;  y  =  0.5;  right  tip. 


Fig.  3.17  Normalized  Mode  II  stress  intensity  factor  vs.  distance 
from  asperity  centerline:  c  =  3a;  \i  =  1.0;  left  tip. 


Fig.  3.18  Normalized  Mode  II  stress  intensity  factor  vs.  distance 
from  asperity  centerline:  c  =  3a;  \i  -  1.0;  right  tip. 


Until  recently,  no  suitable  failure  criterion  could  be  found  for 
the  mixed  mode  fracture.  Instead,  attention  has  been  given  only  to  the 
problem  of  predicting  the  direction  of  crack  extension  when  a  body  with 
cracks  is  loaded.  There  are  two  major  criteria:  maximum  hoop  stress  [58] 
and  minimum  strain  energy  density  [59].  The  maximum  hoop  stress  criterion 
states  that  crack  growth  will  occur  in  a  direction  perpendicular  to  the 
maximum  principal  stress.  On  the  other  hand,  the  minimum  strain -energy- 
density  factor  criterion  postulates  (1)  that  the  initial  crack  growth 
takes  place  in  the  direction  along  which  the  strain-energy-density  factor 
possesses  a  stationary  value,  and  (2)  that  crack  initiation  occurs  when 
the  factor  reaches  a  critical  value.  Comparison  of  the  two  criteria 
by  several  authors  has  shown  that  for  tensile  loading  the  differences 
between  them  are  small.  For  compressive  loading,  however,  not  only 
these  two  do  not  agree  well  but  also  neither  criterion  correlates  well 
with  physically  observed  behavior  [51]. 

When  these  criteria  are  applied  to  the  subsurface  cracks  in  a 
compressive  zone  they  predict  crack  extension  direction  to  be  about  -70° 
(110®  from  x-axis)  at  the  left  tip  and  about  70®  at  the  right  tip 
from  the  direction  parallel  to  the  surface,  implying  that  crack  extension 
occurs  toward  the  surface  at  both  tips.  However,  experimental  results 
by  Jahanmir  [41]  show  that  the  subsurface  cracks  grow  parallel  to  the 
surface  most  of  the  time  before  they  become  loose.  McClintock  [52]  has 
suggested,  in  the  investigation  of  crack  behavior  in  the  rail  head  under 
rolling  condition,  that  cracks  in  a  compression  field  are  most  likely  to 
grow  in  shear.  In  fact,  Forsyth  [60]  has  observed  that  fatigue  cracks 


98 


« 


have  two  growth  regimes.  In  Stage  I,  cracks  formed  on  the  slip  planes 
of  the  persistent  slip  bands  grow  when  they  are  most  closely  aligned 
with  the  maximum  shear  stress  directions. 

The  crack  growth  In  the  shear  direction  has  been  observed  and 
suggested  In  many  cases.  In  the  study  of  ductile  fracture  with  rotation 
Tipnis  and  Cook  [61]  have  shown  that  ductile  fracture  occurs  In  the 
shear  zones  by  either  the  tall  mechanism  when  an  Inclusion  consists  of 
a  hard  core  and  two  soft  tails,  or  the  void  mechanism  when  Inclusion- 
matrix  decohesion  or  brittle  fracture  of  an  Inclusion  occurs,  because 
most  engineering  materials  contain  Inhomogeneities  which  provide  nucleating 
points  for  the  fracture.  While,  Smith  [62]  and  Shleh  [63]  were  able  to 
find  the  compressive  shear  cracks  which  propagate  In  the  slip  bands. 
Further,  Besuner  [64]  has  suggested  that  If  the  cyclic  principal  stresses 
are  primarily  compressive,  then  there  Is  no  reason  for  a  crack  to  align 
Itself  at  right  angles  to  a  compressive  maximum  principal  stress. 

Therefore,  It  Is  necessary  to  investigate  the  variation  of  the  maximum 
shear  stress  with  angle  e,  which  Is  treated  below. 

In  2-d1mens1onal  deformation  field  the  stresses  at  the  crack  tip 
are  expressed  as  [58] 

a  =  — - — [Kt  cos|(3  -  cose)  +  K,,  s1n4{3  cose-  1)] 

2y?¥r  ^  ^  11  ^ 

o.  =  — ^  cos|{K,  cos^  -  y  K,,  sine]  (3.1) 

®  ^  2  I  2  2  II 

1  0 

=  — ■ —  cosy[KT  sine  +  Kt,(3  cose  -1)] 

2^  2  1  II 
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and  the  maximum  shear  stress  t  is  given  by 

max  •' 

^  sin^e)]^^^  (3.2) 

^max  maximum  values  when  the  conditions  of  ^ 

2  2 

3  ®  satisfied.  If  we  denote  0^  for  the  angle  e  which 

satisfies  the  conditions  and  substitute  into  Eq.  3.1,  then  the  stresses 

become  [a„Ja  »  [«Ja]a  »  and  [t^.].  respectively.  Thus,  using  the  Mohr's 
®m  ®  *’m  ®m 

circle  transformation,  the  maximum  of  t  ^  and  the  angle  between 

hiqX 

0  -direction  and  this  maximum  shear  direction  can  be  determined.  When 
m 

the  distance  r  goes  to  zero,  this  direction  ultimately  becomes  the 
direction  of  crack  propagation. 

If  the  above  criterion  is  applied  to  the  results  obtained  in  the 
previous  section  for  the  subsurface  cracks,  it  predicts  the  angle  to  be 
between  -5  and  5  degrees.  These  values  are  very  small  and  therefore  can 
practically  be  assumed  zero,  implying  that  cracks  propagate  in  a  plane 
coincident  with  the  original  cracks  parallel  to  the  surface. 

3.5  Crack  Tip  Plasticity  Consideration 

Whether  the  linear  elastic  fracture  mechanic  (LEFM)  approach  is 
appropriate  depends  on  the  plastic  zone  size  at  the  crack  tip.  If  this 
size  is  too  big  or  comparable  with  such  dimension  as  the  distance  from 
the  crack  tip  to  a  free  surface,  LEFM  is  no  longer  valid  and  plasticity 
plays  a  significant  role.  In  this  case  the  local  stress  and  strain 
history  may  become  important  in  determining  the  fracture  criterion  [65]. 


TOO 


/ 


The  estimate  of  the  plastic  zone  size  for  Mode  II  loading  Is 
given  by 


211  '  k  ^ 


(3.3) 


where  AKjj  Is  the  stress  Intensity  factor  range  In  Mode  II  and  k  Is  the 

shear  yield  strength  of  the  material.  Tables  3.1  and  3.2  show  the 

changes  In  stress  Intensity  factor  and  estimated  plastic  zone  sizes  as 

a  function  of  the  depth  of  crack  location  for  a  small  crack  during  one 

cycle  of  loading.  Values  for  all  parameters  are  the  same  as  given  before 

2 

k  =  245  MPa  =  25  kg/mm  ,  a  =  10  urn.  From  these  two  tables  It  can  be 
seen  that  the  plastic  zone  sizes  estimated  are  comparable  with  the  depth 
of  crack  and  the  stress  Intensity  factor  ranges  calculated  are  close  to 
or  less  than  the  threshold  Intensity  [66].  Of  course,  there  Is  no 
data  In  Mode  II  available.  However,  an  approximate  analysis  treated 
In  Appendix  E  Indicates  that  the  threshold  Intensity  factor  In  Mode  II 
Is  very  close  to  that  In  Mode  I. 

There  are  two  difficulties  In  accepting  the  linear  elastic  fracture 
mechanics  approach  to  the  subsurface  crack  problem.  The  cracks  quite 
near  the  surface  can  have  large  values  of  aKjj,  but  they  cannot  propa¬ 
gate  In  a  brittle  manner  due  to  the  development  of  large  plastic  zones. 

On  the  other  hand,  cracks  far  below  from  the  surface  cannot  grow  since 
the  stress  Intensity  factors  for  these  are  much  smaller  than  the  thres¬ 
hold  although  the  plastic  zones  are  small  and  the  plasticity  restriction 
does  not  apply.  Therefore,  plasticity  should  come  Into  the  analysis 
of  subsurface  crack  behavior  under  the  moving  asperity. 
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Table  3.1.  Stress  Intensity  Factor  Range,  AKjj  (MNtn"^^^),  at 
Left  and  Right  Tips  for  a  Small  Crack  (c  =  l/4a) 


d/a 

p:  0.25 

0.5 

1.0 

Left 

1.05 

1.07 

1.44 

0.3 

Right 

1.28 

1.34 

1.54 

Left 

1.08 

1.11 

1.15 

0.5 

Right 

1.35 

1.40 

1.49 

Left 

0.67 

0.68 

0.70 

Right 

0.89 

0.94 

1.06 

a  =  half  length  of  asperity  contact 
d  =  depth  of  crack  location 
\i  =  coefficient  of  friction 
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Table  3.2.  Ratio  of  Plastic  Zone  Size  to  Depth  of  Crack  at 
Left  and  Right  Tips  for  a  Small  Crack  (c  =  l/4a) 


d/a 

u:  0.25 

0.5 

1.0 

Left 

0.96 

1.00 

1.48 

0.3 

Right 

1.44 

1.56 

2.08 

Left 

0.76 

0.80 

0.88 

0.5 

Right 

1.20 

1.28 

1.48 

Left 

0.12 

0.12 

0.12 

1.0 

Right 

0.20 

0.24 

0.32 

In  summary,  linear  elastic  fracture  mechanics  is  found  to  be  useful 
in  assessing  the  crack  tip  stress  concentration  and  in  determining  the 
crack  trajectory  by  mixed  mode  fracture  criteria.  However,  the  size 
restriction  by  plasticity  consideration  and  the  small  stress  intensity 
factors  calculated  (which  is  less  than  the  threshold)  suggest  that  the 
plastic  fracture  mechanics  approach  is  required. 
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4.  THE  MECHANISM  OF  SUBSURFACE  CRACK  PROPAGATION  - 
PLASTICITY  APPROACH 


4.1  Introduction 

It  was  shown  in  the  previous  chapter  that  linear  elastic  fracture 
mechanics  is  inadequate  in  explaining  the  mechanism  of  subsurface  crack 
propagation.  In  this  chapter  the  mechanism  of  crack  propagation  is 
investigated  through  the  elastic-plastic  finite  element  analysis.  The 
first  part  of  it  deals  with  the  nature  of  plastic  deformation  under 
moving  load  condition  when  there  is  no  crack.  Then,  the  elastic-plastic 
response  of  a  crack  in  the  subsurface  is  studied  and  the  mechanism  of 
crack  propagation  is  discussed. 

4.2  The  Elastic-Plastic  Finite  Element  Model 

In  this  study  AOINA  (Automatic  Dynamic  Incremental  Nonlinear  Analysis, 
a  finite  element  computer  program  for  the  static  and  dynamic  displace¬ 
ment  and  stress  analysis  of  solids,  fluid-structure  systems  and 
structures)  was  used  as  in  the  elastic  analysis  case  in  Chapter  3  [67]. 

The  material  model  used  in  the  study  was  as  follows:  infinitesimal, 
material  nonlinear  only  [68].  The  isothermal  elastic-plastic  incremental 
constitutive  relations  are  shown  in  Appendix  F. 

The  model  used  to  calculate  the  elastic-plastic  response  under  the 
moving  load  condition  was  the  same  as  that  for  the  elastic  case.  No 

dynamic  effect  was  considered  in  the  analysis.  The  material  was  assumed 

-4 

to  be  slightly  work- hardening  (Ej  =  10  E).  The  material  properties 

used  were  as  follows:  isotropic,  slightly  work-hardening. 
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E  =  1.96  X  10^  MPa  =  2  X  10^  kg/nm^,  Ej  =  19.6  MPa  =  2  kgW, 

V  =  0.28,  0^  =  3.^  =  424  MPa  =  43.3  kg/nm^. 

The  finite  element  mesh  used  for  the  study  of  plastic  deformation 
without  a  crack  was  the  same  as  shown  in  Fig.  3.5.  For  the  investi¬ 
gation  of  crack  propagation  only  a  short  crack  (c  =  l/4a)  shown  in 
Fig.  3.4  was  used.  Due  to  the  prohibitively  expensive  computer  cost 
only  a  limited  parameter  study  was  conducted  for  the  case  of  a  =  10  pm, 

Pq  =  4k,  and  y  =  0.25. 

The  problem  was  solved  incrementally  by  moving  the  load  step  by 
step.  For  an  accurate  solution,  the  load  increments  per  step  should  be 
sufficiently  small.  However,  such  a  load  step  requires  a  large  number 
of  calculations  that  make  the  analysis  very  expensive.  Therefore,  larger 
load  steps  were  used  with  iteration  to  obtain  efficient  and  accurate 
solutions. 

The  use  of  iteration  can  introduce  some  difficulties.  The  convergence 
process  may  be  slow,  requiring  a  large  number  of  iterations  that  can  be 
expensive.  Also,  some  iterative  methods  do  not  converge  for  certain 
types  of  problems  or  large  load  increments. 

In  ADINA  program  the  following  procedures  are  implemented  for  the 
iterative  solution  of  the  incremental  equilibrium  equations:  (a)  modified 
Newton  iteration  with/without  Aitken  acceleration  and  a  divergence 
procedure  and  (b)  the  Broyden-Fletcher-Goldfarb-Shanno  (BFGS)  method. 

These  are  described  in  detail  in  Reference  69  by  Bathe  and  Cimento. 

The  modified  Newton  iteration  procedure  is  based  on  the  stiffness  matrix 
at  time  t,  ^K,  while  the  BFGS  method  updates  the  matrix  based  on  the 
iteration  history.  When  divergence  is  detected  for  the  modified  Newton 
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method,  the  divergence  scheme  is  incorporated,  which  uses  smaller  load 
increments.  According  to  Bathe  and  Cimento  iteration  divergence  can 
occur  whenever  the  system  stiffens  during  the  solution  increment. 

For  the  present  study  the  BFGS  method  was  used  because  it  has  better 
convergence  properties. 

4.3  Numerical  Results 

4.3.1  Plastic  Deformation  and  Residual  Stresses 

Fig.  4.1  shows  the  plastically  deformed  region  under  a 
moving  asperity  during  the  first  and  fourth  cycle,  respectively.  Under 
the  moving  load  the  material  just  in  front  of  the  load  gets  plastically 
deformed  and  some  part  of  the  plastic  region  behind  it  becomes  unloaded. 
Due  to  this  internal  unloading  the  shape  of  the  plastically  deformed 
region  under  moving  load  is  quite  different  from  that  under  static 
load  (see  Fig.  4.2).  Furthermore,  the  repeated  loading-unloading  makes 
the  plastic  region  smaller  as  shown  in  Fig.  4.1(b).  According  to  the 
numerical  investigation  of  a  rolling  contact  by  Anand  [70],  a  steady 
state  after  a  few  revolutions  of  a  disk  would  eventually  reach  a  purely 
elastic  state  for  a  strain-hardening  material  whereas  elastic-plastic 
steady  state  condition  could  be  expected  for  elastic-perfectly  plastic 
materials  at  higher  loads.  Therefore,  the  case  shown  in  the  figure  may 
eventually  reach  a  purely  elastic  steady  state  since  the  material  is 
slightly  work-hardening. 

The  residual  stress  is  given  in  Fig.  4.3  as  a  function  of  depth 
from  the  surface.  It  can  be  seen  that  the  variation  in  is  almost 
negligible  after  3  or  4  cycles.  Also,  the  dominant  residual  stress  is 
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Fig.  4. 


1  Plastically  deformed  region  under  a  moving  asperity: 
(a)  during  the  first  cycle;  (b)  during  the  fourth 
cycle. 
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Fig.  4.2  Plastic  zone  under  a  stationary  asperity. 
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compressive  in  a  direction  parallel  to  the  surface.  The  other  components 

of  residual  stresses,  o J"  and  t  are  very  small,  which  should  be  zero 

j  xy 

to  satisfy  equilibrium  [100].  Since  these  two  are  the  only  components  that 
can  affect  crack  opening  and  sliding,  it  can  be  concluded  that  residual 
stresses  do  not  affect  crack  propagation. 

Due  to  the  compressive  residual  stress  the  stress  is  always 
compressive  right  behind  the  asperity  contact  except  very  near  the  surface. 
If  there  is  a  crack  running  perpendicular  to  the  surface,  the  crack  may 
grow  downward  only  a  small  amount  of  less  than  one  quarter  of  the  asperity 
contact  length.  Then,  it  will  cease  to  grow  further  downward  because  of 
the  compressive  stress  that  prevents  the  crack  opening.  Therefore,  it 
may  be  considered  that  the  magnitude  of  the  residual  stress  is  large 
enough  to  make  the  steady  state  stress  compressive,  which  prevents 
cracks  from  running  downward. 

4.3.2  Crack  Tip  Field  Under  Moving  Load 

Proper  crack  tip  modeling  is  quite  important  in  simulation 
of  crack  tip  deformation  and  propagation.  As  discussed  in  the  previous 
chapter  a  triangular  element  has  a  1/r  singularity  when  three  nodes  are 
at  a  crack  tip,  being  allowed  to  separate  during  deformation,  and  the 
side-nodes  are  at  the  midpoints  (see  Fig.  3.1b).  This  type  of  singularity 
representation  is  suitable  for  the  HRR  field  [71-73]  when  the  material 
model  is  elastic-perfectly  plastic. 

In  the  field  of  finite  element  analysis  of  structures  Nagtegaal 
et  al.  [74]  have  found  that  some  frequently  used  elements,  including  the 
three-noded  triangle  and  the  four-noded  isoparametric  element,  are  not 
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suited  for  analysis  in  the  fully  plastic  range.  For  example,  load- 
displacement  curves  obtained  by  finite  element  calculations  generally 
exceed  known  limit  loads,  and  in  some  cases  the  solutions  have  no  limit 
load  at  all.  Also,  deLorenzi  and  Shih  [75]  have  shown  that  the  eight- 
noded  elements  in  general  have  the  same  deficiency.  The  exact  limit 
load  for  a  beam  in  bending  was,  however,  obtained  when  eight-noded 
isoparametric  elements  with  straight  edges  were  used.  Further  study 
by  the  same  authors  suggests  that  the  element  coupled  with  the  appropriate 
near-tip  configuration  possessing  a  1/r  singularity  is  ideally  suited 
for  elastic-plastic  crack  analyses  [76]. 

When  the  plastic  elements  mentioned  above  are  used  in  the  case  of 
tensile  loading,  they  allow  the  modeling  of  crack  tip  blunting  as  the 
nodal  points  separate  during  deformation.  However,  when  a  crack  tip  is 
in  the  combined  compression  and  shear  loading,  all  nodes  at  the  tip 
must  be  constrained  properly.  Without  any  constraint  too  many  crack 
tip  nodes  may  result  in  a  geometrically  incompatible  state.  Results  of 
preliminary  calculations  for  different  crack  tip  modeling  have  suggested 
that  only  one  element  with  1/r  singularity  in  the  direction  of  propagation 
and  several  singular  elements  of  other  type  around  the  element  should 
be  a  good  mesh  design.  This  is  shown  in  Fig.  4.4. 

If  the  crack  tip  is  constrained  in  such  a  way  that  all  nodes  move 

together  during  deformation  (i.e.,  if  all  crack  tip  elements  are  the  type 

of  l//r  singularity),  the  equivalent  strain  at  the  nearest  integration 

points  can  reach  as  high  as  about  200  for  a  crack  located  at  d  =  0.5a. 

(3x3  integration  order  was  used.)  In  Figs.  4.5  and  4.6  the  shear 

strain  y  is  plotted  as  a  function  of  distance  from  the  crack  tip  under 
xy 


Fig.  4.4  Finite  element  mesh  around  a  crack  tip. 
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moving  load  for  both  leading  and  trailing  tips.  The  shear  strains  are 
initially  low  when  sliding  begins  and  reach  very  high  values,  remaining 
nearly  at  these  levels.  It  can  be  seen  that  the  shear  strain  at  the 
left  tip  is  about  twice  as  high  as  the  strain  at  the  right 

tip.  If  the  strains  are  so  large,  then  the  material  at  tips 
cannot  support  such  strains  and  fracture  will  occur  immediately.  Although 
this  model  shows  that  the  strains  at  the  tip  can  reach  infinity,  it 
cannot  adequately  describe  the  details  of  crack  tip  deformation. 

The  plastic  singular  element  can  allow  the  modeling  of  crack  tip 
shear  deformation  when  it  is  used  for  shear  loading.  Fig.  4.7  shows  the 
development  of  plastic  zone  along  with  the  moving  load  step  by  step  when 
this  element  is  used.  At  the  beginning  the  shape  of  the  plastically 
deformed  zone  is  more  or  less  the  same  as  shown  in  Fig.  4.2  for  the  case 
of  no  crack  inside  except  right  around  crack  tips.  Due  to  the  presence 
of  a  crack  the  stress  field  changes  very  much  when  the  load  moves  over 
the  crack.  It  can  also  be  noticed  that  there  are  some  spots  inside  the 
plastic  region  where  unloading  has  taken  place.  With  the  repeated 
loading-unloading  it  is  expected  that  the  overall  plastic  zone  should 
become  smaller  as  in  the  case  of  no  crack  shown  in  Fig.  4.1. 

Examination  of  the  displacements  of  the  nodal  points  on  the  crack 
surfaces  shows  that  the  upper  surface  initially  slides  forward  and  then 
it  slides  backward.  This  is  consistent  with  the  result  by  McClintock  [52] 
who  investigated  the  behavior  of  a  subsurface  crack  in  rolling  contact 
problem.  In  Table  4.1  the  relative  sliding  displacements  of  crack  tip 
nodel  points  are  shown  for  two  different  depths  of  crack  location.  It 
shows  that  the  relative  displacement  increases  with  sliding,  then 
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Fig.  4.7  (Cont'd) 


Table  4.1.  Relative  Crack-Tip  Sliding  Displacement  (um) 
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decreases,  and  finally  changes  the  sign.  Followings  can  also  be 
observed  from  the  table:  (1)  the  crack  closer  to  the  surface  has  large 
values,  and  (2)  the  left  tip  usually  has  large  values  of  relative 
sliding  displacements. 

In  Figs.  4.8  and  4.9  the  crack  tip  shear  strains  are  plotted  along 
the  distance  from  the  crack  tip  for  different  stage  of  the  loading 
position  as  the  asperity  contact  moves  from  left  to  right.  The  distri¬ 
bution  of  shear  strain  increases  until  it  reaches  a  maximum  value  and 
then  decreases.  After  it  attains  a  minimum,  it  increases  again.  At 
very  near  the  tip  the  shear  strain  changes  from  positive  to  negative, 
and  then  to  positive  again.  For  a  crack  located  at  d  =  0.25a,  the 
shear  strain  can  reach  more  than  120^  at  a  point  of  0.00625  ym  away  from 
the  tip. 

4.4  Discussion  on  Crack  Propagation 

The  elastic  solutions  for  cracks  have  a  stress  singularity  of  the 
type  1/y^.  The  singular  terms  of  the  solutions  are  valid  within  a  region 
of  radius  r^  near  a  crack  tip.  Because  of  the  singularity  in  the  stresses, 
the  yield  condition  is  quickly  reached  at  the  crack  tip.  If  the  plastic 
deformation  is  limited  to  a  region  of  radius  rp  and  distor¬ 

tion  of  stresses  on  r^  caused  by  the  plastic  deformation  within  rp  is 
negligible.  In  this  case  the  stress  distribution  in  the  elastic  region 
is  uniquely  specified  by  the  stress  intensity  factors  and  the  crack 
propagation  is  completely  characterized  by  the  critical  stress  intensity 
factor.  However,  the  results  of  the  plastic  zone  estimation  in  Chapter  3 
and  the  elastic-plastic  stress  calculation  in  the  previous  section  clearly 
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Fig.  4.8  Shear  strain  vs.  distance  from  the  left  crack  tip 
when  plastic  elements  are  used,  (a),  (b),  (c), 
(d),  (e),  (f),  (g),  and  (h)  refer  to  the  position 
of  the  left  tip  relative  to  the  asperity  contact 
shown  in  Fig.  4.7. 
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indicate  that  LEFM  is  no  longer  valid  for  the  present  sliding  contact 
problem  and  plastic  fracture  mechanics  must  be  used. 

4.4.1  Fracture  Parameters  beyond  Linear  Elastcity: 

A  Review 

In  recent  years  considerable  effort  has  been  devoted  to 
extend  fracture  mechanics  into  the  elastic-plastic  regime.  These  started 
by  examining  the  stress  distribution  under  fully  plastic  conditions  and 
small-scale  yielding.  Among  several  fracture  parameters  the  J  integral  and 

the  crack  opening  displacement  (COD)  have  attracted  the  attention  of  many 
researchers  since  these  are  viable  candidates  from  the  computation  and 
measurement  point  of  view  [77]. 

For  nonlinear  elastic  material,  it  is  possible  to  express  the 
potential  energy  release  rate  as  the  line  integral,  J,  [78,  79], 
evaluated  around  a  contour,  r,  which  circumscribes  the  crack  tip  in  the 
counter-clock  wise  sense,  starting  from  a  point  on  the  bottom  edge  of  a 
crack: 

f  ^^i 

J  =  J^(W  dy  -  T.-^  ds)  (4.1) 

where 

W  =  the  strain  energy  density  =  o.  .de--j 

J  1 J  1 J  > 

=  the  traction  vector  defined  according  to  the  outward  normal 
n  along  r,  T^.  =  n^, 

Ui  = 


S 


displacement  in  the  direction  of  T^, 
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ds  =  incremental  arc  length  along  r,  and 
x,y  =  rectangular  coordinates. 

The  properties  of  J,  for  nonlinear  elastic  material,  are  (1)  that 
it  is  path  independent,  (2)  that  it  characterizes  the  crack-tip  stress 
and  displacement  fields,  and  (3)  that  it  is  a  measure  of  the  elastic 
energy  release  rate,  G.  For  linear  elastic  material,  J  is  identical  to  G. 

Considerable  work  on  developing  the  J  integral  as  an  analytical  tool 
for  plastic  fracture  mechanics  has  been  performed  by  Begley  and 
Landes  [80,  81]  using  a  compliance  technique.  Following  this  many  investi¬ 
gators  have  been  involved  in  fracture  characterization  by  J.  For  example, 
almost  all  the  papers  in  the  book.  Elastic  Plastic  Fracture,  ASTM  STP  668, 
are  devoted  to  J  integral  approach.  Among  them  some  authors  have  shown 
that  J  has  a  linear  relation  with  COD  [82-84]. 

The  crack-opening  displacement  (COD),  first  proposed  by  Wells  [85] 
in  1961,  is  another  powerful  tool  in  elastic-plastic  fracture  mechanics. 

The  suggestion  of  Wells  is  essentially  equivalent  to  assuming  that  an 
intensity  of  plastic  strain,  critical  for  fracture,  exists  adjacent  to 
the  crack  border  [65],  leading  to  a  near  tip  strain  fracture  criterion  [86]. 
The  COD  concept  has  been  supported  by  several  authors  experimentally 
and  analytically  [87-90].  Analytical  models  as  Pelloux's  alternating 
shear  model  [87]  and  Kuo  and  Liu's  unzipping  model  [90]  indicate  that 
crack  growth  is  approximately  one  half  of  the  crack  tip  opening  displace¬ 
ment  if  rewelding  at  a  crack  tip  does  not  occur.  Depending  on  the  included 
angle  of  the  deformed  tip  crack  growth  can  reach  COD. 


4.4.2  The  Crack  Propagation  Model  In  Wear 

As  reviewed  briefly  the  J  integral  and  the  crack-opening 
displacement  are  the  most  frequently  used  parameters  in  characterization 
of  crack  tip  behavior  in  plastic  fracture  mechanics.  However,  some 
difficulties  are  associated  with  these  parameters.  First,  the  path 
independence  of  J  integral,  which  is  derived  using  deformation  theory  of 
plasticity,  has  not  been  proved  analytically  for  incremental  plasticity. 

With  regard  to  the  characterization  of  the  crack  tip  field,  McClintock  has 
shown  that  the  different  degrees  of  constraint  associated  with  various 
in-plane  configurations  would  prevent  the  satisfactory  development  of  a  one- 
pa  rametertheory  of  plastic  fracture  mechanics  for  full  plasticity  [101]. 
Elastic-plastic  finite  element  calculations  have  shown  that  J  is  strongly 
path-dependent  for  contours  very  close  to  the  crack  tip,  which  is  the  region 
where  heavy  plastic  deformation  has  taken  place  [91,  92],  Therefore,  the 
J  integral  is  not  appropriate  for  the  present  study. 

On  the  other  hand,  the  crack-opening  displacement  is  conceptually 
simple  and  straightforward  regardless  of  its  appropriateness  for  a  one 
parameter  characterization  of  crack  tip  fields,  which  is  still  under 
investigation  by  many  researchers  [77,  93].  The  analytical  treatment, 
however,  is  very  difficult  because  of  the  complexity  for  elastic-plastic 
materials  and  therefore,  finite  element  analyses  have  been  used 
extensively  in  studying  crack  tip  profiles.  In  finite  element  analysis 
it  is  essential  to  use  sophisticated  crack  tip  elements  to  obtain 
satisfactory  modeling  of  the  deformation.  In  addition  to  the  computa¬ 
tional  difficulty,  the  interpretation  of  the  crack-opening  displacement 
in  relation  to  crack  propagation  is  still  to  be  advanced. 

So  far  almost  all  the  crack  tip  deformation  studies  are  limited  to 
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Mode  I  and  Mode  III.  In  sliding  wear  the  subsurface  cracks  are  subjected 
to  Mode  II  loading.  The  cracks  are  likely  to  grow  in  the  direction 
parallel  to  the  surface.  This  has  already  been  proved  in  Chapter  3 
even  though  it  was  based  on  elastic  solutions.  Also  the  experimental 
findings  of  surface  texture  development  due  to  sliding  [94,  95]  strongly 
suggest  that  the  cracks  propagate  along  the  slip  planes  nearly  parallel 
to  the  surface.  Of  course,  there  is  a  possibility  that  crack  growth 
should  follow  the  maximum  principal  stress  or  the  minimum  strain-energy- 
density  factor  criterion  when  a  crack  is  in  tensile  region.  However, 
the  existence  of  a  tensile  component  alone  is  not  sufficient  to  cause 
crack  growth  in  the  direction  that  those  criteria  predict,  unless  the 
strength  of  stress  concentration  is  large  enough.  Later  on  this  pos¬ 
sibility  will  be  discussed  further. 

The  crack-opening  displacement  concept  may  equally  be  applied  to 
the  cracks  in  Mode  II.  In  this  mode  the  relative  sliding  displacement 
occurs  by  means  of  slip  due  to  crack  tip  deformation.  If  the  maximum 
of  this  displacement  is  employed  as  a  crack- tip  sliding  displacement 
(CTSD),  denoted  by  AS,  then  the  crack  growth  length,  AC,  may  be 
expressed  as 


AC  =  AS  -  AC^  (4.2) 

w 

where  AC  is  the  length  of  rewelding.  As  for  aC^  there  is  no  satisfactory 
criterion.  According  to  Kikukawa  et  a1.  [88]  rewelding  is  affected  by 
environmental  conditions  for  the  Mode  I  case.  However,  in  Mode  II  they 
found  that  the  length  of  rewelding  and  therefore,  the  crack  growth  length 
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does  not  seem  to  differ  discernably  between  the  air  and  vacuum  tests. 
Moreover,  the  ratio  of  aC/aS  was  found  to  be  nearly  0.16,  which  is  fairly 
small  in  comparison  with  the  ratio  of  aC/COD  (-x-O.SS)  for  Mode  I  in  their 
study. 

As  implied  by  Kikukawa  et  al.'s  study  a  crack  can  grow  to  the  point 

where  the  fracture  conditions  are  satisfied.  These  conditions  may  be 

addressed  in  terms  of  fracture  strain,  If  the  material  element  at  a 

distance  r^  from  the  crack  tip  attains  the  fracture  strain,  actual  fracture 

will  occur  up  to  this  point.  The  strain  y  along  the  distance  r  from  the 

xy 

crack  tip  can  be  written  as 


where  is  the  strain  intensity  factor,  r^  is  a  characteristic  distance 
used  for  non-dimensional ization,  and  m  is  the  parameter  which  depends  on 
the  strain  hardening  exponent.  For  a  non-hardening  material  in  the  HRR 
field  m  is  unity.  Therefore,  r^  may  be  obtained  by  substituting  for 
Y^y  E'l*  ^-3 


Tf  = 

Fig.  4.10  shows  the  shear  strain  distribution  calculated  for  different 
depths  of  crack  location  as  a  function  of  the  distance  from  the  left  tip. 
From  the  figure  we  may  argue  that  if  y^  is  0.4  or  so,  then  the  cracks  extend 
by  M).004  ym  for  d  =  5  ym  and  'vO.03  ym  for  d  =  2.5  ym. 

When  a  crack  grows  continuously,  it  is  expected  that  the  deformation 


13C 


field  will  change.  Analytic  investigations  of  extending  cracks  are 
limited  to  Mode  III  conditions  [65,  79,  96]  because  of  mathematical 
complexity.  For  Node  I,  only  the  general  functional  form  of  the  near 
tip  singular  strain  field  is  known  for  elastic-perfectly  plastic 
materials  [79].  These  studies  have  shown  that  the  strain  field  ahead 
of  an  extending  crack  is  dominated  by  a  logarithmic  singularity  which 
is  weaker  than  the  1/r  singularity  experienced  at  the  tip  of  a  station¬ 
ary  crack.  Finite  element  solutions  to  extending  crack  problems  include 
simulation  of  crack  advance  by  using  a  nodal  release  technique  [97-99]. 

In  this  method  nodes  are  released  along  the  crack  extension  path  in  a 
predetermined  way,  thereby  propagating  the  crack  the  length  of  one  or 
more  elements.  The  main  objective  of  the  simulation  work  is  to  identify 
possible  macroscopic  parameters  which  characterize  the  deformation  field 
at  the  growing  crack  tip.  This  type  of  fracture  characterization  is 
beyond  the  scope  of  the  present  issue. 

In  summary,  the  propagation  of  subsurface  cracks  in  wear  can  be 
characterized  by  such  fracture  parameters  as  the  crack-tip  sliding 
displacement  and  the  fracture  strain.  For  accurate  numerical  results 
it  is  essential  to  use  specially  designed  crack  tip  elements.  The  iso¬ 
parametric  triangular  elements  with  1/r  singularity  seem  to  be  well 
fitted  to  the  requirements.  Unfortunately,  this  element  is  only  for 
elastic-perfectly  plastic  material. 

As  noticed  the  present  investigation  is  limited  to  the  case  of 
p  =  0.25  due  to  the  prohibitively  expensive  computer  cost.  For  extensive 
study  of  crack  propagation  a  low  cost  method  is  needed.  In  the  study  of 
crack  behavior  in  the  rail  head  under  the  rolling  condition  McClintock  [52] 


has  shown  that  the  boundary  integral  method  is  effective  with  reasonable 
cost.  For  his  problem  the  residual  stress  at  the  shakedown  limit,  super¬ 
imposed  on  a  repeated  rolling  load  small  enough  for  small  scale  yielding 
around  a  horizontal  crack  buried  at  the  point  of  maximum  shakedown  stress, 
was  used  for  the  boundary  conditions.  Here,  we  wish  to  consider  stresses 
above  the  shakedown  limit,  so  the  plastic  zone  extends  to  infinity  in 
the  X  direction.  In  retrospect,  this  could  have  done  with  periodic 
boundary  conditions.  It  should  be  used  in  the  future. 


132 


5.  THE  DELAMINATION  WEAR 


5.1  Introduction 

It  was  shown  in  Chapter  2  that  asperity  deformation,  plowing  by 
wear  particles  and  hard  asperities,  and  adhesion  are  responsible  for 
friction.  These  friction  generating  mechanisms  are  quite  important  in 
understanding  the  delamination  wear  mechanisms.  In  Chapter  3  and  4 
the  subsurface  cracks  in  wear  were  shown  to  be  likely  to  grow  in  a 
ductile  manner.  The  quantitative  model  of  wear  is  described  using  the 
crack  propagation  model  and  the  implication  of  the  mechanisms  of  fric¬ 
tion  and  crack  propagation  is  discussed  in  this  chapter. 

5.2  The  Wear  Model 

Consider  a  subsurface  crack  lying  below  the  surface  as  shown  in 
Fig.  5.1.  An  asperity  is  moving  over  the  surface  from  left  to  right. 
The  wear  rate  will  be  dictated  by  the  crack  propagation  at  both  ends, 

L  and  R.  The  crack  propagation  for  the  i  th  cycle,  AC^. ,  may  be 
expressed  as 


aC^.  =  f(u,  d,  C,  material  properties)  (5.1) 

for  both  ends.  If  N  is  the  total  number  of  asperity  passes  required 
for  removal  of  one  layer,  the  volume  V.|  for  one  crack  of  width  w 
lying  at  a  depth  d  is  obtained  as 

N 

V,  =  W*d»E(AC,  +  ACn  ). 

I  .  L.  R. 


(5.2) 
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crack  under  a  moving  asperity. 
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Therefore,  the  total  volume  V  for  one  layer  may  be  given  by 

N 

V  =  N  •  N  •  w  •  d  •  Z(aC,  +  ACn  )  (5,3) 

c  w  i  “-i  ^ 

where  N  is  the  number  of  cracks  along  the  sliding  direction  and  N  is 
c  w 

the  number  of  wear  sheets  in  the  direction  of  contact  width.  Since 
X  w  is  in  the  order  of  the  contact  width  L^,  the  volume  V  becomes 

N 

V  =  Nc  •  •  d  ‘  z(aCl_  +  aC^  ).  (5.4) 

Let  X  be  the  asperity  contact  spacing,  «.  the  crack  spacing,  D  the 
diameter  of  a  specimen,  and  L  the  contact  length  as  shown  in  Fig.  5.2. 
Using  this  model  the  numbers  and  N  can  be  determined  for  a  specimen 
and  a  slider,  respectively. 

Specimen: 

Since  the  number  of  asperities  per  length  of  contact  is  L/x,  the 
number  N  is  determined  as 


N  -  i_  k 
irD  X 


(5.5) 


where  S  is  the  sliding  distance  required  for  removal  of  one  layer. 
Also,  N  is  given  by 


(5.6) 
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Fig.  5.2  A  model  of  wearing  specimen  and  slider. 
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The  wear  rate  is,  by  assuming  L  =  L,  finally  obtained  as 

y  •  d  •  (AC|^  +  AC|^) 

S  " 


(5.7) 


where  aCj^  and  aCj^  are  the  average  crack  propagation  rates  during  N  cycles. 
Slider: 

The  numbers  N  and  are  given  by 


(5.8) 


and 


(5.9) 


respectively.  Therefore,  the  wear  rate  for  a  slider  is  obtained  as 


V 

S 


(AC^  +  ACp) 


(5.10) 


which  is  the  same  expression  as  given  in  Eq.  5.7  for  a  specimen.  The 
wear  coefficient  K  can  also  be  obtained  by  using  Archard's  equation  as 


K  = 


3H  •  •  d  •  (AC^  +  aC^) 

W  •  X  • 

c 


(5.11) 


where  H  is  the  hardness  of  the  material  and  W  is  the  applied  load. 

According  to  Gupta  and  Cook  [14],  the  number  of  asperity  contacts,  n, 

increases  by  the  0.91  power  and  the  area  of  each  contact  grows  only  by 

the  0.09  power  of  applied  normal  load.  Therefore,  we  may  practically 

assume  the  following  linear  relationship  between  the  real  area  of  contact, 
2 

or  nirr^  ,  and  the  applied  load, 

=  W/H  (5.12) 


where  r^  is  the  average  radius  of  asperity  contact.  On  the  other  hand, 
the  asperity  spacing,  X,  is  proportional  to  L/*^  if  we  assume  a  souare 
array  of  asperity  contacts  and  therefore. 


K 


3d 


ir  •  r  •  £ 

0  c 


(5.13) 


showing  that  K  is  expressed  in  terms  of  the  crack  parameters  (d,  AC's) 
and  asperity  contact  parameters  (r^.  A).  All  the  parameters  involved  in 
deriving  Eq.  5.13  are  summarized  in  Table  5.1. 

At  present  the  validity  of  the  wear  model  described  in  Eq.  5.11  can 
only  be  checked  by  estimating  those  parameters.  Micrographs  of  worn 
specimens  show  that  d  is  of  the  order  of  5  pm  [41].  Estimation  of  A 
from  contact  geometry  was  shown  to  be  of  the  order  of  100  pm.  (Accord¬ 
ing  to  Gupta  [102],  A  is  about  300  pm.)  For  the  values  of  AC's  the  data 
obtained  in  Chapter  4  are  used  and  the  crack  spacing  is  assumed  to  be 

of  the  same  order  as  of  A.  In  the  following  two  examples  the  order  of 
magnitude  estimate  for  K  is  shown  to  be  in  the  ball  park  although  it  is 
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Table  5.1.  List  of  the  Parair«ters  Used  for  Deriving  Eq.  5.13 

X  asperity  contact  spacing 

crack  spacing 

aC  average  crack  growth 

AC^  crack  growth  for  the  i  th  cycle  of  loading 

d  depth  of  crack  location 

n  number  of  asperity  contacts 

r^  average  radius  of  asperity  contact 

w  crack  width 

A  real  area  of  contact 

r 

D  diameter  of  a  specimen 

H  hardness  of  the  material 

L  apparent  contact  length 

apparent  contact  width 

N  total  number  of  asperity  passes  required  for  removal 
of  one  layer 

number  of  cracks  along  the  sliding  direction 
N  number  of  wear  sheets  in  the  direction  of  width 

W 

S  sliding  distance  required  for  removal  of  one  layer 

V  total  volume  of  one  layer  worn 

V.|  volume  of  one  crack  of  width  w  and  thickness  d 
W  applied  normal  load 
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still  premature  to  conclude  anything  further. 

Example  1 : 

The  crack  tip  sliding  diaplacement  for  d  =  5  ym  is  0.0035  ym 
(Table  4.1).  If  aC|^  =  =  0.0035  ym  and  the  apparent  contact  length 

L  =  2  mm  and  both  the  asperity  contact  spacing  x  and  the  crack  spacing 
are  assumed  to  be  100  ym,  then 


I  =  (2  ».  X  10-!l(p.0p7  X  10-‘)  .  ,  ,  ^ 

^  (100  X  10"°)(100  X  10*®) 


and  for  H  =  100  kg/mm  and  W  =  1  kg. 


K  *  3  X  100  X  lO^^x  1.4  X  10*^^  =  4  2  x  10*^ 


Example  2: 

For  d  =  2.5  ym,  by  taking  aCj^  =  ACp  =  0.014  ym  again  from  Table  4.1, 

V  ..  (.3iU°:^)^(2.5  xlO-^)  (0.028X10-^)  ,  ^  3  ^ 

®  (100  X  10*®)  (100  X  100*®) 


and 


K  =  8.2  X  10 


According  to  the  experimental  results  listed  in  Table  2.3,  K  is 


uo 


-2  -4 

between  10  and  10  in  most  cases,  indicatinq  that  at  least  the 
modeling  of  wear  by  the  crack  propagation  study  is  in  the  right  direction. 
The  limitations  on  the  model  are  discussed  below. 

5.3  Discussion  on  Delamination  Wear 

The  delamination  wear  model  described  in  the  previous  section  was 
seen  to  be  able  to  predict  the  wear  behavior  of  materials.  However,  the 
present  study  is  too  limited  to  draw  any  decisive  conclusion.  In  order 
to  generalize  the  model  a  systematic  parametric  study  is  required.  This 
should  include  a  large  amount  of  finite  element  analysis  work  which  is 
prohibitively  expensive. 

Some  factors  such  as  friction  and  fracture  parameters  are  seen  to 
affect  the  delamination  wear  greatly.  High  wear  rates  are  expected  when 
the  coefficient  of  friction  Is  high  and  the  fracture  strain  of  the 
material  worn  Is  low.  Although  our  understanding  of  the  mechanisms  of 
friction  and  crack  propagation  has  been  improved  much,  there  are  still 
certain  limitations  remaining.  In  this  section  the  implications  and 
limitations  of  those  factors  and  some  other  aspects  in  delamination 
wear  are  discussed.  Because  of  the  limited  result  only  qualitative 
discussion  can  be  given.  Both  the  FEM  solutions  for  elastic-plastic 
solids  in  Chapter  4  and  for  elastic  solids  in  Chapter  3  are  used  to 
make  rough  estimates. 

5.3.1  Subsurface  Crack  Propagation 

Eq.  5.7  shows  that  the  wear  rate  is  directly  proportional 
to  the  depth  of  crack  and  the  average  crack  growth.  The  amount  of  crack 
growth  aC  depends  on  the  depth  d.  The  results  in  Chapter  4  indicate 
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that  AC  increases  with  decreasing  d  for  a  given  friction  coefficient 
of  0.25.  For  higher  p  the  same  tendency  is  expected. 

The  reason  for  the  aC  dependence  on  d  may  be  explained  in  terms  of 
shear  strain  concentration  at  crack  tips.  Under  the  moving  load 
condition  it  is  highly  likely  that  the  plastic  deformation  is  much  more 
severe  at  the  crack  tips  near  the  surface  due  to  sliding  between  crack 
surfaces  even  for  low  friction.  However,  the  analysis  of  crack  nuclea- 
tion  [36]  has  shown  that  cracks  cannot  nucleate  at  very  near  the  surface 
because  of  high  triaxial  state  of  compressive  stress. 

For  a  longer  crack  much  larger  shear  concentration  is  expected  due 
to  relatively  easier  sliding  of  crack  surfaces  with  respect  to  each 
other.  Of  course,  this  sliding  may  be  retarded  when  the  shear  traction 
between  the  crack  surfaces  exists.  At  any  rate,  the  crack  propagation 
rate  will  increase  with  increasing  total  crack  length. 

When  the  total  crack  length  reaches  a  certain  size  that  is  critical 
for  unstable  crack  propagation  leading  to  fracture,  the  direction  of 
crack  propagation  will  change  toward  the  surface,  eventually  forming  a 
wear  sheet.  Fracture  criteria  for  combined  mode  of  loading  considered 
in  Chapter  3  may  provide  a  qualitative  explanation  for  this.  Let  us 
consider  a  subsurface  crack  under  a  moving  asperity  shown  in  Fig.  5.1 
again.  If  the  left  tip  L  is  in  the  tensile  zone  behind  the  load, 
then  the  maximum  hoop  stress  and  the  minimum  strain-energy-density 
factor  criteria  predict  the  crack  growth  of  about  -70®  at  the  left  tip. 
When  the  critical  condition  is  reached  with  tensile  component,  the  crack 
will  he  fractured  to  the  surface.  Micrographs  of  wear  sheets  about  to 
del.:^m^nate  show  that  subsurface  cracks  always  reach  the  surface  at  the 
trailing  edge  of  the  crack  [35]. 
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Material  properties  such  as  strain-hardening  property  are  the 
important  factors  affecting  crack  propagation  aC.  Also,  the  anisotropy 
of  the  properties  are  expected  to  affect  aC.  Due  to  repeated  loading 
effects  of  cyclic  hardening,  cyclic  softening  and  the  Bauschinger 
effect  are  also  important.  However,  it  seems  impossible  at  the  present 
time  to  incorporate  all  these  factors  in  the  analysis  because  of  extreme 
complexity.  The  fracture  strain  is  the  simplest  parameter  that  can  be 
used  in  determining  AC. 

5.3.2  Effects  of  Friction  on  Delamination  Wear 

In  order  to  understand  the  delamination  wear  mechanisms  the 
understanding  of  such  friction  mechanisms  as  asperity  deformation,  plowing 
and  adhesion  is  quite  important  for  the  following  reasons: 

(1)  Although  the  overall  coefficient  of  friction  can  reach  a 
large  value,  the  local  surface  traction  at  any  one  of  the  contacts 
cannot  have  a  shear  stress  larger  than  that  predicted  by  plasticity  for 
a  given  type  of  contact.  That  is,  the  shear  stress  at  a  flat  contact 
cannot  exceed  0.4  times  the  normal  stress  if  all  the  normal  load  is 
carried  by  flat  contacts.  Likewise,  the  tangential  stress  exerted  by 
plowing  particles  cannot  be  larger  than  the  normal  stress  since  the 
maximum  value  of  plowing  compont.  t  of  friction  cannot  be  larger  than  1.0. 

(2)  The  deformation  of  asperities  and  the  plowing  action  will 
always  occur,  generating  wear  particles,  since  delamination  wear  keeps 
creating  wear  particles  and  surface  asperities. 

(3)  All  these  three  forms  of  surface  interactions  induce  the 
subsurface  deformation,  crack  nucleation  and  propagation,  eventually 
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causing  delamination  wear.  When  plowing  is  the  primary  cause  of  surface 
traction,  the  surface  will  show  many  grooves  and  remain  rough.  On  the 
other  hand,  adhesion  and  asperity  interaction  will  tend  to  leave  a 
relatively  smooth  surface. 

From  the  computation  point  of  view,  the  above  considerations  can 
raise  an  important  question:  How  can  the  contact  be  modeled  when  the 
coefficient  of  friction  is  much  larger  than  that  possible  in  terms  of 
yield  criteria?  This  is  more  serious  when  the  material  model  is  elastic- 
perfectly  plastic.  There  is  no  satisfactory  answer  to  this  at  present. 
For  a  perfectly  flat  contact,  the  slip  line  solution  indicates  that  the 
normal  stress  decreases  from  (2  +  ii)k  to  (1  +  0.5ii)k  with  increasing 
friction.  This  implies  that  the  real  area  of  contact  should  increase 
twice,  on  which  no  definitive  experimental  evidence  is  available  at  this 
time.  When  the  other  two  components  of  friction  are  present  in  addition 
to  adhesion,  it  is  even  more  difficult. 

Experimentally  it  is  well  known  that  the  wear  behavior  of  materials 
strongly  depends  on  the  coefficient  of  friction  [3,  4].  High  shear 
strain  distribution  is  expected  ahead  of  extending  cracks  when  the 
coefficient  of  friction  is  high.  Therefore,  the  crack  propagation  rate 
is  larger  for  higher  friction  coefficient.  The  exact  relation  between 
friction  and  wear  cannot  be  given  because  of  limited  data.  However, 
the  elastic  solutions  obtained  in  Chapter  3,  although  they  have  already 
been  shown  to  be  inadequate,  may  be  used  to  discuss  the  functional 
relation  qualitatively. 

Tn  Figs.  5.3  and  5.4  the  maximum  Mode  II  stress  intensity  factors 
for  a  o’ven  d  are  plotted  as  a  function  of  friction  coefficient  y. 
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Fig- 


5.4  Normalized  Kjj  as  a  function  of  y 
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It  can  be  seen  that  the  relation  between  Kjj  and  y  is  approximately  linear 
for  a  given  d.  If  we  assume  that  Paris'  crack  growth  law  is  still  valid 
in  this  case,  then  the  crack  growth  AC  may  be  expressed  as 

AC -v-  (av)^  (5.14) 

where  a  is  a  constant.  Since  the  wear  volume  is  determined  by  the 
amount  of  the  total  aC  accumulated,  wear  rates  are  a  power  function  of 
friction.  The  exponent  b  is  a  constant  which  varies  from  2  to  4  depend¬ 
ing  on  material  for  aK  values  in  stage  II.  This  should  give  a  rough 
idea  about  a  functional  relation  between  friction  and  wear,  although  the 
foregoing  discussion  is  very  approximate  and  may  not  be  valid  for 
subsurface  crack  propagation  in  wear. 

When  aK  is  less  than  the  value  in  stage  II  or  the  threshold  value, 
normally  this  is  the  case  when  the  applied  load  is  small,  Eq.  5.14  is  not 
applicable.  For  the  case  we  may  assume  the  following  functional  form: 

f=f(rp)  (5.15) 

where  rp  is  the  pi  Stic  zone  size  given  by  Eq.  3.3.  To  determine  the 
function  f  the  strain  distribution  in  the  pi  Stic  zone  and  the  fracture 
strain  should  be  known.  However,  it  does  not  seem  to  be  possible  now. 

If  we  assume  the  following  linear  form  as  a  first  approximation. 


(5.16) 


then  the  exponent  b  in  Eq.  5.14  is  2  and  Eq.  5.15  becomes 


dN  2ii  ^  k  ^ 


(5.17) 


where  a  is  a  constant.  Since  &Kjj  =  'for  a  given  d,  Eq.  5.17  can  be 

written  as 


dC  a  /Pox2  f. 


(5.18) 


and  int'igrating  both  sides  after  separation  yields  to 


C  =  exp 


(5.19) 


where  and  3  are  constants.  Therefore,  the  wear  volume  V,  using 
Eq,  5.4,  is  obtained  as 

V  =  •  d  .  2Cp  exp  (-j^)  -N].  (5.20) 

All  the  parameters  used  for  deriving  Eq.  5.20  are  summarized  in  Table  5.2. 
To  use  this  expression  to  estimate  the  wear  volume,  all  the  constants 
should  be  determined  first  and  more  fundamentally,  the  exact  form  of  f  in 
Eq.  5.15  should  be  known, 

5.4  ConcU: 'ons 

(1)  The  wear  equation  in  the  delamination  theory  is  promising  in 
predicting  the  wear  of  materials. 

(2)  The  subsurface  crack  propagation  in  wear  can  be  determined 
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Table  5.2.  List  of  the  Parameters  Used  for  Deriving  Eq.  5.20 

a,  3  proportionality  constants 
AKjj  Mode  II  stress  intensity  factor  range 

a  proportionality  constant 

b  exponent 

d  depth  of  crack  location 

k  shear  yield  stress  of  the  material 
Pq  maximum  stress  value  at  the  center  of  asperity  contact 
Cq  constant  of  integration 

apparent  contact  width 

N  total  number  of  asperity  passes  required  for  removal 

of  one  layer 

number  of  cracks  along  the  sliding  direction 
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by  such  plastic  fracture  parameters  as  crack-tip  sliding 
displacement  and  fracture  strain. 

(3)  The  wear  rate  of  material  is  expected  to  be  a  power  function 
of  the  coefficient  of  friction  when  the  test  conditions  are 


identical . 
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6.  SUMMARY  AND  CONCLUSIONS 

The  experimental  and  theoretical  investigation  reported  in  this 
study  verified  the  hypothesis  of  the  delamination  theory  of  wear  which 
states  that  wear  occurs  as  a  result  of  force  transmission,  plastic 
deformation,  crack  nucleation,  crack  propagation,  and  wear  sheet  formation. 
The  theoretical  anlyais  indicated  that  the  wear  rate  of  metals  can  be 
predicted  by  friction  and  crack  propagation  information. 

Studies  were  carried  out  on  the  friction  mechanisms,  showing  that 
the  coefficient  of  friction  between  the  sliding  surfaces  is  due  to  the 
combined  effects  of  asperity  deformation,  plowing  by  wear  particles  and 
hard  surface  asperities,  and  adhesion  between  the  flat  surfaces.  It  was 
found  that  the  frictional  behavior  of  materials  is  history  dependent 
and  can  successfully  be  represented  in  the  "Friction  Space"  diagram  shown 
in  Fig.  2.18.  Further,  the  canpatibility  of  sliding  surfaces  was  found  to 
be  governed  as  much  by  mechanical  properties  of  materials  as  by  the  adhe¬ 
sion  between  the  two  surfaces.  This  investigation  indicated  that  to 
incorporate  all  the  effects  of  friction  components  into  modeling  of 
sliding  contact  for  deformation  study  further  research  is  required, 
particularly  for  the  case  of  complicated  contact  geometry. 

Studies  on  the  crack  propagation  in  wear  indicated  that  a  subsurface 
crack  grows  in  the  direction  parallel  to  the  surface  until  a  critical 
state  is  reached  at  the  trailing  edge  of  the  crack,  turning  the  direction 
toward  the  surface.  Considerations  on  the  plastic  zone  at  the  crack  tip 
and  the  threshold  stress  intensity  factors  suggested  that  the  linear 
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elastic  fracture  mechanics  approach  is  inadequate  for  subsurface  cracks 
in  sliding  wear.  As  expected  from  the  result  of  Johnson  and  Jefferis  [103] 
a  semi -infinite  solid  under  the  cyclic  moving  load  was  found  not  to  reach 
shakedown  for  p^  =  4k  and  u  =  0.25  if  the  material  model  is  elastic- 
perfectly  plastic.  Such  fracture  parameters  as  crack-tip  sliding  displace¬ 
ment  and  fracture  strain  were  found  to  be  useful  for  characterization  of 
crack  tip  field.  For  a  complete  characterization  of  crack  propagation 
further  research  is  required  in  the  field  of  fracture  mechanics  as  well  as 
in  wear. 

Finally,  the  delamination  wear  equation  given  as  Eq.  5.11,  which  was 
derived  using  the  crack  propagation  model,  was  found  to  be  promising  in 
predicting  the  wear  of  materials.  Due  to  limited  data  in  the  crack  propa¬ 
gation  study  for  the  elastic-plastic  solids,  anything  decisive  such  as 
the  exact  relation  between  friction  and  wear  cannot  be  determined. 

However,  the  rough  estimate  by  the  calculated  stress  intensity  factors 
indicated  that  there  is  a  power  functional  relation  between  friction  and 
wear  although  the  stress  intensity  factor  approach  was  shown  to  be  inade¬ 
quate.  When  the  applied  load  is  small,  the  plastic  zone  size  estimated 
by  the  stress  intensity  factor  may  be  useful  in  deriving  the  crack 
growth,  and  therefore,  the  wear  volume,  which  was  shown  in  Eq.  5.15 
through  Eq.  5.20.  For  this  purpose,  further  research  should  be  done  to 
determine  the  relation  between  crack  growth  rate  and  plastic  zone  size. 
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APPENDIX  A:  SLIP-LINE  FIELD  SOLUTION  FOR  SLIDING  CONTACT 

Fig.  A.l  gives  a  possible  slip-line  field  for  the  asperity  contact 
between  sliding  surfaces.  The  interface  OA  between  asperities  and  the 
stress  free  surface  OD  are  both  assumed  to  be  straight  with  their 
directions  defined  by  the  angles  e  and  0'  measured  from  the  sliding 
direction.  From  the  figure  it  can  be  noticed  that  the  slip-line  ABCD 
is  a  6-line.  Using  the  Hencky  relations  the  stresses  along  the  slip¬ 
line  can  be  obtained. 

At  D,  ())  =  0'  +  ^  and  p  =  k.  (A.l) 

Along  AB,  41  =  a  -0  +  |-  and 

Pab  =  k  (1  +  J  +  2a  -  20  -  20')  (A. 2) 

Isolating  the  junction  along  ABO  as  shown  in  Fig.  A. 2  we  can  find  the 
resultant  forces  as 

-Fy  =  L  =  (AB)  p  cos  (0-a)  -  (AB)k  sin  (0-a) 

+  (0B)k  cos  (6-a)  -  (0B)p  sin  (0-a)  (A. 3) 

-Fx  =  F  =  (AB)p  sin  (6-a)  +  (AB)k  cos  (0-a) 

+  (0B)k  sin  (0-a)  +  (0B)p  cos  (6-a). 

Using  the  geometric  relations,  L  and  F  are  expressed  as 

L  =  (OA)  (p  cose  +  k  sin  (2a-0))  (A. 4) 

F  =  (OA)  (p  sine  +  k  cos  (2a-0)). 
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Therefore,  the  coefficient  of  friction  is 

c  (1+  |-  +  2a-  2e  -  26')  sine  +  cos  (2a-9) 

v=T  = - i -  (A. 5) 

(1+  I  +  2a-  2e  -  26*)  cose  +  sin  (2a-e). 

When  e,  e',  and  a  are  small  (e,  e',  o-+0) 

M  =  ~ —  =  0.39. 

1  .f 

For  several  values  of  a,  e,  and  e',  the  coefficient  of  friction  is  plotted 
in  Fig.  2.12. 

If  the  junction  does  not  weld  along  OA  the  shear  stress  along  OA 
will  be  much  less  than  k.  The  interfacial  shear  stress  is  related  to 
the  angle  a  as  T=k  cos  2a.  When  there  is  no  shear  stress  along  OA,  i.e., 
a=TT/4,  the  junction  will  slide  along  OA  until  the  junction  can  deform 
under  the  influence  of  the  normal  load  alone.  During  this  sliding  the 
coefficient  of  friction  can  be  obtained  by  substituting  a  =  7r/4,  into 
Eg.  (A. 5)  as 

(1  +  n-2e-  2e')  sine  +  cos  (|-  -e) 

VI  =  - - - 

(1  +  Tr-2e  -  2e')  cose  +  sin  -e) 
or 

M  -  tan  e 

which  is  the  same  expression  as  that  derived  from  the  roughness  theory 
of  friction. 

The  slip-line  field  solution  derived  above  is  a  general  solution. 

Wh^ '  the  sliding  occurs,  the  slip-line  field  should  satisfy  the  kinematic 
(  ’■‘:ion,  which  corresponds  to  the  case  of  the  slip-line  AB  being 

oarallel  to  the  sliding  direction.  Therefore,  o  is  equal  to  e. 
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APPENDIX  B:  PLOWING  COMPONENT  OF  FRICTION 

The  plowing  force  P  can  be  expressed  as  P  =  pA^  where  p  is  the 

flow  pressure  and  A^  is  the  cross-sectional  area  of  the  grooved  track. 

Sin  et  al.  [B.l]  have  shown  that  the  expression  for  the  plowing 

component  of  friction  Hp  given  by  Goddard  and  Wilman  [B.2]  is  found 

to  agree  well  with  the  experimental  results.  In  this  appendix  the 

expression  is  derived  following  their  treatment. 

Let  us  consider  a  spherical  particle  of  radius  r  loaded  by  a 

2 

normal  force  L.  Then  the  real  area  of  contact,  irr^^,  is  obtained  by 
=  L  (B,„ 

where  r^  is  the  radius  of  the  periphery  of  the  indentation  and  H  is 
the  hardness  of  the  material.  When  a  tangential  force  is  applied,  just 
sufficient  to  cause  very  slow  sliding,  a  groove  is  plowed  out  by  the 
particle,  and  the  indentation  must  become  deeper  to  support  the  same 
load,  because  the  contact  interface  is  then  only  on  the  front  half  of 
the  particle.  If  the  pile-up  above  the  initial  metal  surface  in  front 
of  the  moving  particle  is  not  too  large,  the  horizontal  projected  area 
of  contact  is  virtually  semicircular,  the  radius  w/2  being  such  that: 


where  p  is  the  normal  pressure  that  is  smaller  than  H  in  this  case. 
Therefore,  the  plowing  force  can  be  calculated  using  this  value  of  w, 
and  taking  into  account  that  the  effective  flow  pressure  Pg  applicable 
here  may  be  somewhat  larger  than  p.  From  the  geometry  of  the  plowed 
groove  the  cross-sectional  area  Ag  is  given  by 
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A, '  =1"'’  ?F  -  r  (“-3) 

Since  =  PgAg,  and  using  Eqs.  (B.2)  and  (B.3), 


which  becomes  the  expression  given  as  Eq.  (2.5)  when  -*■  p. 
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APPENDIX  C;  STRESS  INTENSITY  FACTORS  BY  FEM 

Elements  with  the  mid-side  nodes  placed  at  the  quarter  point  posi¬ 
tion  have  been  employed  in  Refs.  [C.l]  and  [C.2]  to  obtain  the  mode  I 
stress  intensity  factor,  Kj,  for  elastic  crack  problems.  It  was  done 
by  extrapolating  a  displacement  parameter  [C.3],  associated  with  the 
crack  opening  profile,  to  the  crack  tip.  However,  it  was  shown  by 
Shih,  deLorenzi,  and  German  [C.4]  that  to  have  consistency  between  the 
interpolating  function  and  Williams'  eigenfunction  expansion  [C.5] 
along  the  edge  ABC  in  Fig.  C.l,  the  stress  intensity  factor  Kj,  must 
be  related  to  the  nodal  displacements  u^,  Ug,  and  Uq  with  Up^=0  in  most 
cases  by 


v  -  -3u^  +  Aug  -  Ug 

•^I  ■  TiT~  ^ 

where  G  is  the  shear  modulus, k  =  3-4v  for  plane  strain  and  k=  (3-v)/ 
(1+v)  for  plane  stress.  For  mixed  mode  loading,  the  stress  intensity 
factors  can  be  obtained  by  relating  the  stresses  given  by  FEM  to  those 
by  eigenfunction  expansion.  This  is  treated  here. 

The  geometry  of  an  8-noded  plane  isoparametric  element  is  mapped 
into  the  normalized  square  space  (t,  s),  (-If^t^l,  -IfS^l)  through  the 
transformations  (Fig.  C.2), 

8 


y  = 


8 

1=1  ^  ^i 


(C.2) 


1  Isoparametric  quarter  point  elements 
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where  h^.  are  interpolation  functions  [C.6]  and  x^. ,  are  the  nodal 
coordinates  of  the  element.  For  the  element  shown  in  Fig.  C.2a  we 
have 


^1  ^  ^2  ""  ^3  ^4  ^  ^6  ^5  “  4’  ^8  ^  f 


substituting  Eq.  (C.3)  in  Eq.  (C.2)  by  using  interpolation  functions, 
and  collecting  terms  we  get 


x  =  J(l+t)2  (C.4) 

y  =  I  (1  +  s) (1  +  t)^ 

Any  point  P  on  the  radial  line  R,  is  at  a  distance  r  from  the  crack 

tip, 

r  =  (x2  +  y2)V2  =  b  ^^^^2  |.(2aj^  ^  (l+s)^]^'"^ 

or 

1/2  2  1/4 

1  +  t  =  /F/  {(|)  [(f^)  +  (1  +  s)2]  }  (C.5) 

The  Jacobian  of  transformation  and  its  inverse  are,  respectively, 
given  by 


■  3X  3^' 
3t  3t 

f(i+t)  |(i+s)(in) 

3X  ^ 

.as  as. 

0 

and  (C.6) 


170 


[0] 


-1 


'  2  1 

_4 

1+s 

a  1+t 

"  a 

0 

8 

1 

F 

(1+t)^ 

The  strains  are  given  by 

1 


■<§  >  <0-7) 

au  ^  av 
ay  ax 

The  derivatives  of  displacements  u  and  v  with  respect  to  x  and  y  above 
are  obtained  by  multiplying  by  their  derivatives  with  respect 

to  t  and  s,  i.e. , 


ly.) 

au 

ax  1 

at 

>=  [jj-'' 

au  1 

au 

ay) 

las) 

and 


(C.8) 


av 

'  av 

ax 

-  [jj"'' 

at 

av 

av 

ay 

as 

The  stress  intensity  factors,  Kj  and  Kjj,  are  normally  defined 
along  the  line  9  =  0  (s  =  -1).  Using  s  *  -1  for  simplicity,  the 
strains  along  this  line  are  given  by 
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^  =  1+t 


0 

0 


2t-l 


^  0 


l  0 


2t+l 


2(-t-2)  0 

b 


2t-l  2(-t-2)  2t+l  4(l-t) 


4(1^1 

b 


a 


•4(l-t)  -4t 
b  a 


0  0 


^  0 


8 

b 


^  0 


lu.} 


(C.9) 


where  {u^}  is  the  displacement  vector.  The  stresses  in  this  element  are 
related  to  the  element  strains  using 

T=Ce  (C.IO) 

where  C  is  the  elasticity  matrix.  For  plane  strain  conditions,  the 
stresses  are 


a 

r- 

o 

> 

> 

1 

' 

XX 

1  /« 

E 

V  1-v  0  , 

XX 

Vy> 

yyy 

>  (l+v)(l-2v) 

V 

0  0 

>, 

or 
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From  Eqs.  (C.5)  and  (C.11)  it  can  be  seen  that  stresses  can  be  written 
In  the  form: 


—  (E.  +  F. .) 


(C.12) 


where  E^ •  and  F. .  are  the  factors  to  be  determined  from  the  result  of 
the  finite  element  analysis.  In  two-dimensional  deformation  field,  the 
stress  distribution  in  the  vicinity  of  a  crack  tip  has  been  shown  to 
take  the  form  [C.7]  : 

K  K 

=  — —  COS  |-  [1  -  sin  j  sin  |^] — ^  sin  |C2+cos  4  cos  |^] 

/2irr  /2TTr 


COS  4  [1+sin  4  sin  |^]  + 


/2Trr 


cos  TT  sin  Tj-  cos  ^ 


.  •'l 


xy 


sin  j  cos  cos  2 


30 


1I_ 


cos  f  [1  -  sin  ^  sin  |^]  (C.13) 


where  Kj  is  the  opening  mode  of  stress  intensity  factor  and  is  defined 
as 


Kj  =  Aim  o  (r,  e=0)  •  y^iir  (C.14) 

and  Kjj  is  the  stress  intensity  factor  for  mode  II  and  is  also  obtained 
by 


K,,  =  Aim  T  (r,  e=0)  •  y^irr  (C.15) 

Therefore,  the  stress  intensity  factors,  from  Eqs.  (C.ll),  (C.14)  and 


(C.15),  can  be  expressed  as 
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K,  =  .^  F. 


and 


yy 


^  fxy 


(C.16) 


Further  manipulation  using  Eq.  (C.5)  with  s  =  -1  leads  to  the  following 
expressions  for  Kj  and  Kjp 

(C.17) 

K„  ■  /Bi  (t^ _,) 


where  q.  .  .  and  o,-,-  7  a>"e  the  stresses  at  nodal  number  4  and  7,  respective- 
ly  (see  Fig.  C.2a).  The  similar  expressions  to  those  in  Eq.  (C.17)  have 
already  been  used  by  Harrop  [C.8]. 
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APPENDIX  D:  FRICTION  BOUNDARY  CONDITION  BETWEEN  CRACK  SURFACES 


When  a  crack  is  located  in  the  compression  region  the  upper  and 
lower  crack  surfaces  will  contact  each  other.  If  there  appears  friction 
between  these  surfaces  the  relative  sliding  motion  with  respect  to  each 
other  will  be  retarded.  For  the  rectangular  coordinate  system  in 
which  the  x-axis  is  parallel  to  the  crack,  the  boundary  conditions  of 
closed  crack  surfaces  may  be  expressed  by  the  following  relations: 


xyt 

=  "'xyu. 

xy 

=  ±MOy 

“  % 

(D.l) 


wnere  the  subscripts  i  and  u  indicate  the  lower  and  the  upper  surface* 

respectively.  The  sign  on  the  right  hand  of  the  third  equation 

should  be  chosen  properly  according  to  the  signs  of  t  and  a  . 

xy  y 

In  the  finite  element  program  based  on  the  displacement  method, 
the  stresses  are  expressed  by  a  linear  combination  of  the  nodal 
displacements  for  the  linear  elastic  case.  Therefore,  the  stress 
boundary  conditions  in  Eq.  (D.l)  can  be  transformed  into  linear 
functions  of  nodal  displacements.  By  choosing  the  appropriate  relation 
the  frictional  boundary  condition  problem  can  be  solved  by  FEM. 
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APPENDIX  E:  THRESHOLD  STRESS  INTENSITY  RANGE  IN  MODE  II 


For  cracks  in  Mode  II  to  propagate,  aKjj  should  be  greater  than 
ti.5  threshold  range  in  Mode  II.  However,  there  is  no  experimental 
data  or  theoretical  work  on  threshold  in  Mode  II  available  presently. 
A  very  approximate  expression  for  the  threshold  stress  intensity 
factor  in  Mode  II  is  derived  by  formally  following  the  recent  work  in 
Mode  I  by  Yu  and  Yan  [E.l]. 

In  monotonic  loading  case  the  plastic  zone  size  at  a  crack  tip, 
b),  can  be  given  by 


1  ^ 


(E.l) 


where  aKjj  is  the  stress  intensity  factor  range  and  is  the  shear 
yield  stress.  Also,  the  strain  in  the  plastic  zone  'i  is  assumed  to 
follow  a  modified  HRR  field  solution  of  the  form 


*  /  01  \ 
^  -  ^y 


l/(l+n) 


(E.2) 


where  y  is  the  shear  yield  strain,  x  the  distance  from  the  crack  tip, 
n  the  strain-hardening  exponent,  and  p  the  factor  added  to  prevent  the 
plastic  strain  from  reaching  infinity  at  x=0  [E.2].  Thus,  the  average 
plastic  strain  over  a  length  i  is  obtained  by 


'avg  i 


=  -r  r  ^  rJ 


y'  n  ''t 


] 

(E.3) 


For  a  threshold  (at  the  initiation  of  crack  propagation),  the 
process  zone  t  and  also  the  factor  p  may  reach  the  critical  root 
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root  radius  or  the  Burgers  vector  |Bl  since  the  basic  deformation 

process  of  metals  is  slip.  By  assuming  reaches  the  fracture 

flvg 

strain  Y|:  and  by  substituting  Eg.  (E.1)  into  Eq.  (E.3),  we  obtain 

ak2 

T,  =  r  - pS-)  [jn/d+n)  .  ,]  (E.4) 

“min 

and  therefore  from  this  can  be  solved  as 


AK 


Gy 

nth  “  V 


(E.5) 


where  G  is  the  shear  modulus.  By  taking  n=1  as  a  first  approximation, 
we  get 


“llth  ‘  ’-2  (E-6) 

If  we  assume  y^  =  [2(l+v)//5]e^,  then  with  G  =  E/2(l+v)  and  v  =0.3 
the  threshold  becomes 


^•^Ilth  "  °'^^®f’^’'Pmin 


(E.7) 


which  shows  =  0.7 

References : 

E.1.  C.  Yu  and  M.  Yan,  "A  Calculation  of  the  Threshold  Stress  Intensity 
Range  for  Fatigue  Crack  Propagation  in  Metals,"  Fatigue  Eng.  Mat. 
Str. ,  Vol.  3,  1980,  pp.  189-192. 


E.2.  S.  d.  Antolovich  and  A.  Saxena,  "A  Model  for  Fatigue  Crack 
Propagation,"  Eng.  Fract.  Mech.,  Vol.  7,  1975,  pp.  649-652. 


i 


ill 


179 


APPENDIX  F;  ELASTIC- PLASTIC  INCREMENTAL  CONSTITUTIVE  REUTIONS  USED 
IN  FEM  ANALYSIS  [F.l] 

In  elastic-plastic  analysis  the  material  behavior  is  described 
using  three  properties  in  addition  to  the  elastic  stress-strain 
relations,  C^:  (1)  a  yield  condition,  (2)  a  flow  rule,  and  (3)  a 
hardening  rule. 

The  initial  and  subsequent  yield  condition  for  isothermal  kinematic 
or  isotropic  hardening  can  be  written  as 

^F(V,  V  ^<)  =  0  (F.l) 


where  for  isotropic  hardening  the  strain  hardening  parameter  depends 
on  the  plastic  deformation,  and  for  kinematic  hardening  is  a  constant. 
Restricting  the  analysis  to  associated  flow  rules,  the  function  F  in 
Eq.  (F.l)  is  the  plastic  potential  function  to  which  the  normality 
rule  is  applicable,  i.e.. 


'ij 


‘a  ’‘f 


3% 


IJ 


or  in  matrix  form 


(F.2) 


P 

e 


(F.3) 


P  t 

where  e^j  is  the  plastic  incremental  strains  and  x  is  a  scalar  to  be 

determined.  Since  during  plastic  deformation  ^F=0,  we  also  have 


d^F  = 


3^F 

3  a.. 


Oi-  . 


3^F 

3%P. 

IJ 


=  0 


or  in  matrix  form 


(F.4) 


£  e 


(F.5) 
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where  o--  is  the  incremental  stresses.  The  stress  increments  are 

*  J 

calculated  from 


a  = 


(F.6) 


Using  Eqs.  (F.3)  and  (F.6)  to  eliminate  £  and  a,  we  obtain 


3.  £  £ 


‘x= 


(F.7) 


Substituting  in  Eq.  (F.6)  from  Eqns.  (F.3)  and  (F.7)  for  e^*^  and  ^x. 
the  elastic-plastic  material  law  at  time  t  becomes 

(F.8) 


trEP  ^ 
a  =  C  e 


where 


tnT  t„  .  t  T  pE V 
£  £  +  £  £  £ 


F.9) 


In  ADINA  the  von  Mises  and  the  Drucker-Prager  yield  condit  ons  have 
been  implemented.  When  the  von  Mises  yield  condition  is  used  the 
loading  surface  for  isothermal  classical  isotropic  or  kinematic  hardening 
is  given  by 


tp(  a,  e  ,  k)  -  ^(  s^.j  -  ot^j)(  s^j 


t  \ 
-  “ij^ 


(F.IO) 


where  ^s^j  is  the  deviatoric  stresses  and  is  a  tensor  denoting 
the  translation  of  the  yield  surface 

In  the  case  of  isotropic  hardening,  which  includes  perfectly  plastic 
conditions,  we  have 


t  2 
a.. 


/3 


(F.ll) 
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where  is  the  yield  stress  at  time  t  and  is  a  function  of  the  plastic 

t  P 

work  per  unit  volume  W  , 


t  P 

e,.. 


(F.12) 


Evaluating  and  we  obrain  in  two-dimensional  analysis 


(F.13) 


and 


V 


=  'h 


ft 

[  Oil 


°11  °12 


‘^33^ 


where 

=  2/3 


(F.U) 


(F.15) 


We  can  now  calculate  the  material  property  matrix  in  Eq.  (F.9).  In 
Table  F.l  the  stress-strain  relations  for  two-dimensional  plain  strain 
analysis  are  given.  In  the  case  of  linear  hardening,  if  the  strain¬ 
hardening  modulus  Ej  (tangential  modulus)  is  available,  can  be 
calculated  as 


1/^H  =  3/2  (1/E.p  -  1/E) 


(F.16) 
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